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PREFACE 


It is widely felt that when the equations of mathematical 
physics are written out in full Cartesian form the structural 
simplicity of the formulae is often hidden by the mechanical 
labour of writing out every term explicitly. Attempts have 
been made to reduce this labour by one form or another of 
vector algebra; but it has always seemed to me that this 
method both introduces new difficulties and is insufficiently 
general. Thus the product of two vectors, in vector lan¬ 
guage, means one of two things, either the scalar or the 
vector product, and it is not physically obvious why just 
these functions of the vectors should arise and no others. 

The use of tensor notation, with the summation con¬ 
vention, carries out as great a simplification of the writing 
as does vector notation. The notation has actually at¬ 
tracted attention owing to its applications in the theory of 
relativity, but for ordinary purposes two great abbrevia¬ 
tions may be made. We use rectangular Cartesian axes; 
the result is that the distinction between covariant and 
contravariant vectors disappears, and with it the terms 
arising from curvature of the surfaces of reference. The 
formidable character of most of the formulae of the theory 
of relativity is absent from the formulae of tensors referred 
to Cartesian axes. The tensor method is a necessity for re¬ 
lativity; for applications in dynamics, electricity, elas¬ 
ticity, and hydrodynamics it is a great convenience. 

It is found that the scalar and vector products are not 
the only functions of two vectors that arise, though the 
theory provides reasons why they are important in many 
applications. There is also a symmetrical product, which 
ordinary vector notation is completely unable to express. 
In tensor notation it arises naturally as a symmetrical 
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tensor of the second order. The system of moments and 
products of inertia of a rigid body constitutes such a 
tensor; so do the stress components and the strain com¬ 
ponents in an elastic solid. 

The present method, like vector notation, is of use prin¬ 
cipally in proving general theorems. In concrete applica¬ 
tions there is usually some asymmetry about the coor¬ 
dinates that makes it necessary to abandon the tensor 
form at some stage in the work. It has been said that 
vector equations are like a pocket map, and it has been re¬ 
plied that a pocket map has to be taken out of the pocket 
and unfolded before it is of any use. The same applies to the 
tensor method, and for the same reason; but it has the 
great advantage that it is not a new notation, but a concise 
way of writing the ordinary notation, so that the unfolding 
can be carried out more conveniently when occasion arises* 

What is usually called Statics is treated in Chapter v, 
after Dynamics. I consider this to be the proper order, be¬ 
cause Statics is a special case of Dynamics, and many of its 
formulae have physical significance for reasons explained 
in Dynamics. The customary reversal of the order is due, 
I believe, to the fact that an introduction to mechanics has 
to be given at schools before the students have received 
any training in calculus; but this need not influence 
students working for a university examination* 

It should perhaps be stated that the object of this work 
is to illustrate the use of tensor methods; it does not claim 
to give a complete theory of all the subjects touched, re¬ 
ference for which must be made to the standard text-books. 

I must express my gratitude to Mr M. H* A* Newman, 
Miss L* M. Swain, Dr 3, Goldstein, and Dr Bertha Swiries 
for assistance at various stages in the work, and to the 
staff of the University Press for their care in the printing. 

HAROLD JEFFREYS 

September 1931 


NOTE 


Since this book was written most of the material in it 
has been incorporated in Chapters 2 and 3 of Methods of 
Mathematical Physics, by my wife, Bertha Swiries Jeffreys, 
and myself. The present reprint has been made partly 
because of a continuing demand for a treatment of Car¬ 
tesian tensors by themselves; partly because some results, 
notably on the thermodynamics of an clastic solid and the 
circulation in viscous flow, are not given in textbooks of 
the special subjects. 

HAROLD JEFFREYS 


January , 1952 


NOTE ON THE SEVENTH IMPRESSION 


In the present reprint a few small changes have been made. 
In Ex. 3, p. 15, the components arc given explicitly; these 
have been found useful in developing the st rain energy for 
an elastic sphere under rotation. Omissions have been 
corrected on pp. 8 i and 82. As these do not affect the 
later work they appear to have given readers no trouble. 


December, 19GS 


HAROLD JEFFREYS 








CONTENTS 


Chap. I. 

Cartesian Tensors 

page 1 

n. 

Geometrical Applications 

16 

in. 

Particle Dynamics 

24 

IV. 

Dynamics of Rigid Bodies 

2!) 

V. 

Equivalence of Systems of Forces 

42 

VI. 

Continuous Systems 

48 

TO 

Isotropie Tensors 

66 

VIII. 

Elasticity 

71 

IX. 

Hydrodynamics 

83 

Index 

93 


CHAPTER t 


CARTESIAN TENSORS 


If we have two sets of rectangular axes (Ox, 0y % 0z) t 
(0x\ 0y\ Oz') at the same origin, the coordinates of a 
point P with respect to the second set are given in terms 
of the coordinates with respect to the first set by the 
equations 

x r = l x z + why + rc-iZ’ 

y* = l^x + m^y + n%z • (I). 

z f — l$x + m^y + 7hz i 


The quantities (£*, m t > th, ... a are the cosines of the 
angles between the various axes; thus b is the cosine of the 
angle between the axes Ox ' and Ox, n 2 is the cosine of the 
angle between Oy* and 0z f and so on. It follows that the 
coordinates (x, y\ z) can be expressed in terms of (x' 3 y% z’) 
by the relations 


x = bx* -f by* 4- l % z* 

y = + wi 2 z' - 

z = + n z y' + n z z f , 


(2). 


We can shorten the writing of (I) and (2) considerably 
by a change of notation. Instead of {x, y , z) let us write 
(x lf x 2J Xg) 3 and instead of (z\ y\ z') write ( x t \ x 2 \ x 3 '), We 
can now say that the coordinates with respect to the first 
set of axes are x i} where i may be 1, 2, or 3; and those with 
respect to the second set are x /, where j may bo 1, 2, or 3. 
Then in (i) each coordinate x/ is expressed as the sum of 
three terms depending on the three x i . Each x % is associated 
with the cosine of the angle between the direction of that 
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Xf increasing and that of x/ increasing. Let us denote this 
cosine by a#. Then we have, for all values of j, 

x/ = a u Xi + a 2i x 2 + OzjXt 

= S a fj x t (3). 

<-1.3,3 

Conversely (2) can be written 

x t = S a u x/ (4), 

J-l, 2, 3 

the a u having the same value as in (3), for the same values 
of i and j, because it is in both cases the cosine of the angle 
between the directions of and x/ increasing. 

In mathematical physics we often have to deal with sets 
of three quantities in relation to a set of axes, of the general 
form Ui {that is, 2 ^, %, %), and such that in relation to a 
different set of axes the corresponding quantities are 
(%', u/ y u/) t which satisfy the relations 

u/ = E a H Ut (5) 

<«1, 2. 3 

and u t = E a a uf (6). 

jf“l 2 3 v 1 

Such sets of three quantities are called tensors of the first 
order , or vectors . The individual %, may be called the 
components of the tensors. 

Clearly if we multiply all of the u t and u/ by the same 
quantity m we get 

mu/ E {mu t ) (7), 

i- 1. 2 t 3 

so that mu t is another tensor of the first order* 

Again, if we have two tensors of the first order, u t and 
v it we shall have 

“/ + »/ = s «# i («t + V t ) (8), 

i—1, 2, 3 ' ' 

so that u t + v t is a tensor of the first order. 

We notice that each of the equations (3) to (8) is really 
a set of three equations; where the suffix i or j appears on 
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the left it is to be given in turn all the values 1,2, 3, and the 
resulting equation is asserted in each case. In each such 
equation the right side is the sum of three terms, obtained 
by giving j or i the values 1, 2, 3 in turn and adding. Wher¬ 
ever such a summation occurs a suffix is repeated in the 
expression for the general term; where there is a summa¬ 
tion for all values of j the general term, such as a if u t % con¬ 
tains^ twice . We make it a regular convention that, unless 
the contrary is stated, whenever a suffix is repeated it is 
to be given all possible values and that the terms are to be 
added for all. Thus we write (5) as simply 

u/ = a if Ui (9), 

the summation sign being automatically understood by 
our convention. Then (9) really means three equations, 
with three terms on the right of each, but we can by means 
of our conventions express all of the twelve terms com¬ 
pactly by the single equation (9). 

There are single quantities, such as mass and distance, 
that are the same for all sets of axes. These are called 
ten sors of zero order , or scalars . 

Consider now two tensors of the first order, u t and v k . 
(When we write w a tensor u/' we mean of course a tensor 
of the first order whose components are u lf This is 

another piece of shorthand.) Suppose each component of 
the one multiplied by each component of the other; then 
we obtain a set of nine quantities expressed by UiV k > where 
each of i and k is independently given all the values 1, 2, 3. 
The components of u t , with respect to the other set of 
axes are u/ t v/ say; and 

u/v/ (a a u t ) (a kl v k ) 

= ^ii a kl^iV k ( 10 ), 

The suffixes % and k are repeated on the right. Thus (10) 
represents nine equations, each with nine terms on the 
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right. Each term on the right is the product of two factors, 
one of the form a ti a ki , depending only on the orientation of 
the axes, and the other of the form representing the 
products of the components referred to the original axes. 
In this way the various u/v% can be obtained in terms of 
the original UiV k , But products of two vectors are far from 
being the only quantities satisfying this rule* In general a 
set of nine quantities w ik referred to a set of axes, and 
transformed to another set by the rule 

Wj\ = (I!)* 

is called a tensor of the secovid order. 

We may go on similarly to construct and define tensors 
of the third, fourth, and higher orders. Thus a set of quan¬ 
tities that transforms like XiX k z m x v ... is called a tensor of 
order n r where ti is the number of factors in this product. 

When we say that a certain set of quantities is a tensor 
of any order n, we mean that we have ways of specifying 
its components with respect to any set of axes, and that the 
components with regard to any two different sots of axes 
are related according to the rule appropriate to tensors of 
that order, and in particular to the products of the coor¬ 
dinates with n factors* For instance, if we say that u t is a 
tensor of order 1, we are not simply defining u/ as meaning 
a u Ui. We are supposing both that u/ has a meaning, such 
as a displacement or a velocity, with reference to the axes 
of xf 9 and that the value of each component is equal to 
a u u t . Thus the statement that any set of quantities is a 
tensor is not a mere convention, but a statement capable 
of test and therefore needing proof. In (7) and (8), for 
example, our data are that % and u/ are the components 
of a vector with regard to two different sets of axes* We 
prove that the sets of quantities obtained by multiplying 
both by the same quantity are related according to the 
vector rule; and therefore the products are vectors. 
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If we interchange j and l in (11), we get 

w,/ = a a a kS w ik (12)* 

But on the right i and k are “dummy suffixes”; that is, 
they are to be given all possible values and the results 
added. It is unimportant which of them we call i and winch 
fc; we may therefore interchange them and get 


w ti * = a kl a u w ki = a tj a kl w ki (13), 


Thus w ki transforms according to the same rule as w ik and 
therefore is another tensor of the second order. The im¬ 
portance of tills is that if we know the set of quantities 
arranged 


Mi 

W,2 

Ms 

( *% 

Ml 

Ms 

'Ml 

M* 

Ms. 

s second order, 

/Mi 

Ml 

Mi 

1 Ms 

We 

Ml 

' Mil 

Ms 

Was 


(14) 


( 15 ) 


is another tensor of the second order. Therefore the sets 
(w ik + w ki ) and (iv ik — w ki ) are tensors of the second order. 
The first of these has the property that it is unaltered by 
interchanging i and k f and is therefore called a symmetrical 
tensor. The second has all its components reversed in sign 
when i and k are interchanged, and is called an antisym - 
metrical tensor* Clearly in an anti symmetrical tensor the 
“leading diagonal” components, i,e. those with i and k 
equal, are ail zero. Also, since 

Wik ~ i + I — w m) (Ifr)* 


we can consider any tensor of the second order as the sum 
of symmetrical and antisym metrical parts. 

The gradient of a scalar is a vector. For if 17 is a scalar, 
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its gradient is 3 Ujdx, or dUjdx/ according to the set of 
axes. But 


3 U__ 3 z, dU = SU 
dx/ ~ ax/ dx t a,i dx { 


( 17 ). 


so that the gradients transform according to the vector 
rule. Similarly the gradient of a vector is a tensor of 
order 2. For if u, and «/ are the components of a vector 
with respect to two sets of axes, 


3 uj 

d$i 


dx k du/ 
dx t ' dx k 


a*, ^ («««,) 


du ( 


(18), 


so that the rule of transformation is as in (11). 

Since Xi is a vector, it follows that dx t jdxj. is a tensor of 
the second order. But dx i jdx k is unity ii i = k and zero if 
i ^ k. Hence the set of quantities S (t , such that 



constitutes a tensor of the second order. We can prove this 
directly; for if we apply (11), S j( ' in the new system of 
coordinates should be given by 

Bit = OfiO’kl^tk ( 1 ®)* 

The suffix k has to take all values 1, 2, 3. But if k =£ i, 
S Jfc is 0, and the corresponding term is zero. If k = i, 
= 1, and the result of the summation with regard to k is 

Bji = OffjOti (20)* 

But the a it are the direction cosines of the axis of x/ with 
regard to the x t , and the a n are those of x t ' with regard to 
Xi . Hence a tl a it is the cosine of the angle between */ and 
Xi, and is equal to 1 if the axes are identical and to 0 if they 
are perpendicular. It follows that the result of the trails - 
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formation is that S„' = 1 if j = l, and S,,' = 0 if j & l. It 
follows that the set of quantities 

fl 0 OX 

0 10 ( 21 ) 

\0 0 1 / 

is transformed into itself by the rule (11) and therefore is 
a tensor of the second order. 

If u f is a vector and we form the product S we have 
a tensor of order 3. But now put m = k and add for all 
values of k. Since S lft = 0 except for k = i, the only term 
different from zero is that for k = i, and this is u ( . Hence 

= Ut ( 22 ). 

This operation therefore replaces the suffix k by i. The 
tensor S ik can therefore be called the substitution tensor. 

In the tensor w (k let us put k = i, and in accordance 
with our convention add for all values of i. Then the 
corresponding quantity w^' is got by putting l—j and 
summing; but 

= v>u ( 23 ). 

Thus w it transforms into itself and therefore is a scalar. 

This operation of putting two suffixes in a tensor equal 
and adding accordingly is known as contraction. In general 
it gives a new tensor, whose order is less by 2 than that of 
the original tensor* If for instance we contract the tensor 
u { v kt we obtain 

UiVf - u x v % + UiVz + (24) t 

w‘hich is the scalar product of u ( and v k . 

Similarly the tensor u i v k yields the symmetrical and 
antisymmetrical tensors {u$i k -f u k v { ) and (u i v k ^u k v i ). We 
may call these the symmetrical&nd aiUisymmetrical products 
of u L and v k . 
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The tensor du^dxi gives similarly, on contraction, a scalar 

(25), 


du t _ £% 3 dut + 3% 
dxi ~ dx l dx 2 dx 3 


which is known as the divergence of u ( ; while it gives also 
symmetrical and antisymmetrical tensors 
S u k , du t , du k du t 
dx t + dx k 1 ' 

The former has important applications, especially in the 
theory of elasticity and hydrodynamics; the latter is 
known as the curl or rotation of u { . The vanishing of the 
curl Is the condition that u { may be the gradient of a scalar. 

All the above considerations can be extended to any 
number of dimensions. In n dimensions a tensor of order r 
has n T components, A tensor of order 2, in particular, has 
n ,a components. If it is antisymmetrical s the n diagonal 
components are zero, and the others are equal and opposite 
in pairs. Hence an antisymmetrical tensor of order 2 has 
(n — I) independent components. If ft — 2, 3, 4,,.. in 

turn, this number is 0,1, 3, 6,.... It happens that in three 
dimensions the number of numerically independent com¬ 
ponents of an antisymmetrical tensor of the second order is 
equal to the number of components of a vector. Actually it 
can be proved that with any vector we can associate an anti¬ 
symmetrical tensor of the second order, and conversely. 
This is not true in any number of dimensions other than 3, 
Since the aa are the direction cosines with respect to the x t 
of three perpendicular lines, they are connected by six re¬ 
lations 


a n 2 H- a** 4 — 1 


&& 4 + 

fl m 2 4- a^ 


= 1 , 

+ ffgs&aa. 4 a 32®33 = ^ ! 

fl^Ou + a $i a zi 4 - = 0 

All 013 4* fl^Lflat + 031^32 — 


(20), 


(27). 
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We notice that the second and third of (27) both contain 
(flu, a 2t , a 31 ). We may therefor© solve them for the ratios 
of these quantities. Thus 






023032 ~ a '£i a Z2 a 33 a U ~ G 13^32 ®1 Z a n ~ 0 23^12 


= h 


<*8). 


say. Substituting in the first of (26) we get 

1 — Aj {flu (fl^flgg %%) 4 fl^i ( 0330*2 %%) 

4" 031 {013022 — 023032)} 


k 


Also 


fl n <% 0 31 

012 022 032 

013 023 033 


(29), 


** {(® 23 «a- + (*33*12 - o'-ia^aa)® + {*13*22- *22*12)“} = 1 

(30). 

Rut we have a general identity 

(o* -f 6 s + c a ) (a'* + & #t + c**) - (aa‘ + hb* + cc')* 

= (6c' - cky + (ca f - ac*f + (o6 # - bay (31). 

Hence 

[(013" 4 023^ 4 033 2 ) (012 2 4 022" + 032^) 

— (012013 + 02102B 4 032033)“] = 1 (32), 

Rut on account of the second and third of (26) and the first 
of (27) the expression in brackets is unity, and therefore 

Jfe = ± 1 (33). 

For any given transformation the determinant in (29) is 
therefore equal to ± I. Evidently its sign is reversed if we 
interchange any two of the suffixes j, for this interchanges 
two rows of the determinant; so that the sign is a matter of 
the numbering of the axes. If we start with a rigid frame 
attached to the axes x it and rotate it continuously till it 
is attached to the axes %/, all the a a vary continuously and 
therefore the determinant cannot change from -f- I to — 1 
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or from — 1 to 4- L If then x x goes to x/, x % to x 2 , and x z 
to x zt the determinant ts initially 


and therefore 


1 0 0 
0 1 0 
0 0 1 


1 


Jfc = — I 


(34), 


(35), 


and the determinant formed by the a ti is always + 1. 

If we have a frame of axes (x x x 2 x a ) we can turn it by a 
continuous movement so as to bring x x along the old x zt 
x 2 along the old x a , and along the old x x * In this case we 
have 




ii 

$ 

x % - 

= X 

s> *s' - 

(36), 

and 








a n z 

= 0, 

%1 — If G 3t — 

0; 

0 12 ! 

= 0 

'» «S* = 0> <*32 = 1 

(37). 



a i3 = 

&23 

= 0, 

<*33 = U 

The 

determinant of the c 

hi is 

therefore 





0 

1 

0 

= 1 

(38), 




0 

0 

1 






1 

0 

0 




as before. Any rotation of the axes that does not alter the 
cyclic interchange of suffixes 12 3 12 3 1-.. therefore 
leaves the determinant equal to unity, and therefore so 
long as we always use right-handed or always left-handed 
axes the determinant of the a fi is + 1. 

With this restriction 

fljl — ^ ^32^33 t ^21 = ^13^32 ^33^12 t (30) 

~ ^^12 “ ^ 13^22 J 

and therefore every direction cosine is equal to its first 
minor in the determinant. 
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Those relations are of course identical with those ex¬ 
pressed in the usual notation of solid geometry by 

(40). 

Now suppose that u t is a vector, and consider the set of 
quantities 

Wik = / 0 % - (41). 

( - 0 u t j 

\ 0 / 

Apply (11) to this, taking i to be the number of the row and 
k that of the column. We see that % enters as w n and as 
— Its coefficient in w n ' is therefore a 2i <hi — 
and in all 

w n ~ ((hfthi ~ a zj a w) % + u % 

+ (®i — a 2j<hi) % (42). 

This is obviously zero if j = L If J ^ l and if the other axis 
perpendicular to x/ and x*' is x n * , and j hi jin is a cyclic 
order, the quantities in brackets are equal to [a lni a £n> a Bn ). 
This is true if l immediately succeeds j in the order. If l 
precedes j by one place the signs are reversed. Hence if 
j = I and l = 2, or j = 2 and l ^ 3, or if j = 3 and l = 1, 

w n f = a in u t = u n * (43), 

and in the alternative case 

w n * = - u n f (44). 

Thus w n ' = / 0 u z -V\ (45), 

( - < 0 <) 

\ < -«t' 0 / 

and is of the same form as (41). Thus with any vector we 
can associate an antisymmetrical tensor of the second order. 
Conversely with any antisymmetrical tensor of the second 
order we can associate a vector. 


















<h i 

<hl 

<kl 

<hi 

<*!■ 

®2n 

O-Zn 
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We can proceed alternatively by considering the set of 
quantities defined by the condition that if any two 
of i, k , m are equal the corresponding component is 0; if 
i, k, m are all unequal and in cyclic order, the component 
is + 1; if the order is not cyclic, the component is — 1. 
Let us see whether this is a tensor of the third order. If so* 
we should have 

n — a rt a Al&mn*ikm 

= Q'ifQ'Zl<hn + Q'ZjChl&ln + a 3i^u^2n 

— — &3i (46), 

Now if, for instance, j = l , the right side is clearly zero and 
€nn * = 0, If j , n are all unequal, the expression is 


(47), 


which is equal to 1 if jin are in cyclic order and to — 1 if 
not. Hence the set of quantities e tkm is transformed into 
itself by the rule for transforming tensors of order 3, and 
therefore constitutes a tensor of order 3, This is called the 
alternating tensor. 

Now consider the product € fkm u pi where u p is a vector. 
This is a tensor of the fourth order. If w© contract it by 
putting p = m and summing we get a second order tensor 
w ik — e ikm u m * If i = 1 and k = 2, the only value of m that 
makes e 12m different from zero is 3, and then * ikm = + L Hence 

w n = u 3 (48), 

If i = 2 and k= 1, m is 3; but 213 is the reverse of cyclic 
order and e m = — I, Hence 

w n = “ % t 49 )* 

Similarly we find that the elements of w ik are 
/ 0 — IfA 

( - «3 0 ) 

\ — «! 0 / 


(50), 
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so that the antisymmetrical tensor associated with a vector 
can actually be obtained from it by multiplying by 
and contracting* 

Again, suppose that we are given a tensor of the second 
order w# and that we form a vector u m by multiplying by 
€ ikm and contracting twice. We have, if m = 3, 

= €123^12 4 *213*^21 = ^1E ^21- 

Thus n m = £tkm W tJb (®^)* 

If w ik is symmetrical, this evidently gives zero. If it is 
antisymmetrical the components of u m are numerically 
twice those of w a . 

On account of the intimate relation between the vector 
and the antisymmetrical tensor we shall habitually denote 
the tensor w fk of (41) by u ikJ so that 

= ^22 = — 0 ; a — u lf — w*; j 

Uft = — ^2=“%, 33 — ^2 * 

It will always be seen at once whether the vector or the 
tensor is intended, since the former has one and the latter 
two suffixes. 

If we have any three vectors u it v it w t , and consider the 
scalar € ikm u i v k w mt we see that 

€ ikm UiV k W m = u t v 2 w 3 + u i v 3 w 1 + U 3 V t tV 2 


% 

*» 

% 


V a 


»1 


w a 


so that we have a concise way of writing the determinant 
formed by the components of three vectors. If any two of 
the vectors are parallel this scalar vanishes. 

In associating a vector with an antisymmetrical tensor 
of order 2 a sign convention clearly arises. Wo make the 
positive signs in (41) He one place to the right of the leading 
















u CARTESIAN TENSORS 

diagonal. If then we have two vectors u t and v k , their antU 
symmetrical product is u t Vk — and in the associated 
vector we give the positive sign to w ik when k follows i in 
the cyclic order. Hence the components of this vector are 
taken to be 

{%t?a - u z v t - - n 2 v i). 

This vector is perpendicular to both the original vectors; for 
% - u z v^ + u 2 (%*?* - 4- %{%% - = 0, 

v l (^^3 - *bv 2 ) + % - %^a) + % (% v 2 — u t v t) = °* 

We call it the vector product of % and and can save 
writing by denoting it by [u, v] m . 

Similarly with the antiflymmetrical tensor „ * — we 

associate a vector so as to leave the sign unaltered when 
k follows i in the cyclic order. Thus the components are 

/ dt fa _ 3w» 3% _ 3wg du^ _ 3uj\ ( 

\dx 2 dz 2 5 3a?3 dx k 9 dx v dxj ' 

This is often called “curl u 

The Tensor e ika e mM , Since this tensor is the product of 
two third order tensors, once summed, it is a tensor of 
the fourth order, i t k, m, p being arbitrarily assignable. 
Evidently if { = k or m — p, the corresponding component 
is zero. 

If i = m, the contribution from any value of s is zero 
unless also k = p, and then 

^iks = e mp» = i I> 
and the component is + L 

If i = p, then no value of s gives a contribution unless 
Jfc «s m . Then one of € iki and e mpa is 4 - 1 and the other — 1, 
and the component is — 1* Hence the components of the 
tensor are as follows. 
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If i m t k = p, the component is + I, unless i also = k, 
i = p t k = m, the component is — I, unless i also = k t 
i = k or m = p t the component is 0. 

These results apply also to the tensor 


and therefore 


= Sf m S*j} SfijSini 


(55), 


EXAMPLES 

1. If ti jf u/i u k are the components of a vector with regard to 
three sets of axes, prove that the values of u k * are the same as 
would be obtained by transforming first from u t to u/ arid then 
from u/ to u k 

2. Provo that 1=3 0* 

3. Show that 

( ikm = &a{hk&em ” &ck + 

4. Provo that € in t mkM - 2B im ; - 6, 

5. Provo that *tkM*ms* = «**<*» * - 


6. Provo that if u (J v k * w m are vectors, 
[W [u f v]\ m - u m (v ( w<) - 

[ v * W ]m “ 


7. If 


prove that 


A (») - 


, 


t% 

u 3l 


A (u) «* *imn 


*Wc u il u jm u kn ^fmn A (w), 

6A (w) = *||fc*jffln u U u /m u fcn* 

8, Use Ex. 7 to prove the rule for the multiplication of deter¬ 
minants 


where 


'fmn'imn ^ («) ^ i v ) = Ci («) A ( 1! ) = 6i (“ u )> 
(uv)iv - *hi v lp- 















C H A P T E B II 


GEOMETRICAL APPLICATIONS 

The displacement from any point to any other obviously 
constitutes a vector* The distance between the points is a 
scalar* If x it are the coordinates of the points and r the 
distance between them, 

r % — tit ~ ^) a (*)> 

the square on the right indicating the scalar product of the 
vector into itself. Also the quantities (y t — xj/r constitute 
a vector. 

If we take a fixed point and consider points given by 

= tt* + %iT (^}> 

where r is a variable scalar and the i* are constants such 
that 

V - 1 (3), 

tit - a,)* - r* (4), 

so that r is the distance of a* ( from a t . If we take another 
point yi such that 

y t = a t + hs (5), 

tit - «*)* =* ( 6 )> 

tit - **)* - t* - r)> (7), 

and therefore the distances between a u x if and y, are such 
that the sum of two of them is equal to the third. Thus the 
points are on a straight line; and (2) gives the equations 
of the line in terms of the parameter r. The l t are the 
direction cosines of the line. 

If we take two lines through a, given by 

Xi = a t + h r 
Vl = &i + 


( 8 ), 
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the distance between % and y* is given by 

tii ~ x i ) 2 = ( m i - rl i)* 

= s 2 + r 3 — 2rsl i m i (10). 

Rut this quantity is also equal to s 2 + r 2 — 2 rs cos 9, where 
6 is the angle between the lines. Hence the angle between 
two intersecting lines is given by 

cos 6 liTHi (11). 

If two lines have the same direction cosines they are said 
to be parallel. If two lines do not intersect we can take a 
line through any point on one of them parallel to the other; 
then this line is inclined to the first at an angle given by (11). 
We can then use (11) to determine a unique quantity 
associated with any two lines, which we may call their 
inclination, whether they intersect or not. 

If we have a line given by (2) and y t is a point outside it, 
the line joining and y t subtends a right angle at x t if 

tit - <*i ) 2 = tit - «<)“ + tii ~ x i ) 2 

= rH f * 4 - tit - *< - rl i) 1 

= tit - «<) a “ 2rl i tii - a d + * r *V ( 12 )> 

and therefore 

r^lt tit - a i) ( i3 )- 

This gives the projection of the displacement y t — a t on 
the line. The foot of the perpendicular is 

4- l^r = a* 4 hl k tit ~ &*)• 

Evidently r in (13) will be the same for all points such 
that is constant. Hence 

hVi = * (14) 

represents a plane perpendicular to the line. 

If we take two intersecting lines given by (8) and (0), we 
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can find the equation of the plan© containing them as 
follows. If this plane is 

(15), 

this equation must be satisfied by x { and y, for all values of 
r and 8. Hence 

n i®i — P ( 16 ), 

- 0 (17), 

= 0 (18), 

and from (15) and (16), 

n t (z ( - as,) = 0 (19), 

Then (17), (18), (19) are three homogeneous equations in 
the and can be consistent only if 

e flbn (*i ~ *<) « 0 (20), 

This is the equation of the required plane. Also the n t are 
proportional to 

^ikm^k^'m ~ (4^®9 * 4 ^ 2 * 4 ®H " 4 ^ 3 > 4^2 — ( 21 ), 

But nf = 1 (22), 

(4™a - k™*) 2 + (4™i - h^) 2 + (h m s - 4»h) a 
= (4 3 + 4 2 + V) (^h 2 + + *V) - (4«h + 4^ + 4 ^ 3 ) 2 

= 1 — cos 2 0 

= sin 2 0 (23). 

Thus sin fln ( =± e,**4m m (24). 

The ambiguity in sign corresponds to a general one in 
specifying the parameter r of a point on a line* If r in (8) 
is taken negative, we get a point on the Line on the opposite 
side of a t from those given by positive values of r* But if 
we reverse both r and the l { we still keep If = 1 , and we 
still have the same point* We may take either direction 
along a line to be that of r increasing; if wo reverse the 
direction the signs of ail the 4 are reversed for the same 
point. 
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For any point on (8), 

^ikm x k^m ~ € (km ( a k + 4 r ) 4i = *ikm a k^m (25). 
This is a constant vector for all values of r, and may there¬ 
fore be considered as a property of the line* We denote it 
by 4\ Then we have six properties of the line given by 
4, 44 These are coordinates of the line* They are connected 
by two relations, 

4 a - 1 (26), 

44 = tikm 4a*4t ™ 6 (27). 

The If have a geometrical interpretation* Thus if we con¬ 
sider the plane 

Ss4’^4 = 4' (28), 

this plane passes through the line. Also if x 2 = a 2 , x z ~ 

(28) is satisfied for all values of x lt and therefore if we take 
a line through a t parallel to the x x axis, (28) represents the 
plane through (8) and this line* Two such planes determine 
the line, and therefore the 4 and If together determine the 
line. 

If we have two non-intersecting lines given by 

Xi = a t 4- rl t (29), 

Vi = ^ + sm t (30), 

the line y t = + rl t (31) 

passes through jB* and is parallel to (29), The plane including 
(30) and (31) is, by (20), 

*tkm (2( - Pt) l k m m = 0 (32), 

This therefore represents a plane through (30) parallel to 

(29) , The plane through (29) parallel to (30) is 

*ikm (*< - «<) 4™m = 0 (33). 

The distance between these planes is the projection of the 
line joining any two points on them upon a line parpen- 
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dicular to both. If a line perpendicular to both has direc¬ 
tion cosines n it the shortest distance d between the lines 
is therefore given by 

d sin 9 — (ft — a ( ) n £ sin 9 

= ± (ft — a*) e ikm l k m m 
— ± UikmPi^k^m ~ € ikm a ih m m) (34)* 

But € ikm ^rn m - m k f ; € (kin a f l k = l m ' (35), 

and therefore 

dsin 9 — ± (— l t m/ — ftt f 4')t 

so that, apart from the ambiguity in sign, 

d sin 9 = l t m/ -j- mft' (36)* 

Thus the shortest distance is directly expressible in terms 
of the coordinates of the two lines. 

Now consider two intersecting lines 

x t = a £ + rlii yi=* a £ + (37)* 

The area of the triangle formed by a if x 4 , y t is sin 0. 

The projections of these points on the plane x t = 0 are 
(0, a 2 , £ 3 ), (0, ar 3 ), (0, y 3 ) and form a triangle whose 

area is 


1 

<h 

«* 

= h 

1 


a 3 

1 1 




0 

hr 

V 

I 

y* 

ft 


0 

?n 2 s 

IRsS 


= \rs (4m* - 4 m*) 

— %rst sin $ (38). 

Thus the projections of a triangle, and therefore of any 
plane area, on the coordinate planes are in the ratios of the 
direction cosines of the normal to the planes. A plane area 
can therefore be treated as a vector whose components are 
proportional to the direction cosines of the normal. 
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If we have a line given by 

x £ = at + rlf (39), 

and ft is a point not on the line, let us suppose the point ft 
turned through an angle 6 about the line. The foot of the 
normal from ft to the line is given by 

r - 4 (ft - «<) (40), 

and therefore the displacement from the foot of the normal 
to ft is equal to 

Pt — fat + 44 (ft — %)} = 4 2 (ft — 0*) — 44 (ft — a*) 

— 4 (4 (ft “ ®t) “■ 4 (ft — Gjt)} 

(41)* 

The magnitude of this displacement, p, is given by 

j> 2 -(ft-^) a -{4(ft-^)F 
— 4 2 (Pi — — 44 (ft — (ft ~ «*) 

= i (4 (Pi — & t) — 4 (ft — a*)} 2 (42), 

the | being needed because in the double summation each 
pair of values of the suffixes would occur twice* 

The plane through ft and the line is 

n i (%i “ ®i) = 0, 

subject to nJi = 0, 

n i (Pt “ Gt) = 0, 

and is therefore 

€ ikm (%i — &i) 4 (Pm — ®m) = 0 (43), 

while the ra, are proportional to € ikm l k (ft, - a m ) and there¬ 
fore equal to ± * ikm l k (ft - a m )jp. 

If now we turn ft through an angle 9 , it receives a dis¬ 
placement p (l — cos 9) along the normal to (39), and a 
displacement p sin 9 along the perpendicular to the plane 
(43), If it goes to y i} we have therefore 
Yi - ft « — (1 - cos 9) 4 (4 (ft - a £ ) - 4 (ft - a k )} 

± am 0 € ikm l k (ft - a m ) (44), 
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If 8 is a small angle and we neglect 8\ the displacement is 
simply 

yi — Pi — ± (Pm ~ a m) (46), 

and in particular 

Yt-Pi = ± ihO (Pa (Pa - <h)} (46), 

and from the additive form of tins equation we see that the 
displacement is the sum of those given by separate small 
rotations l { 8 about axes through a t parallel to the co¬ 
ordinate axes. Conversely, displacements due to small 
rotations about axes through a point can be added vectori- 
ally as if all were applied to the system in its original 
position, and give the same total displacement as if they 
were compounded into a single rotation about an axis by 
the vector rule. We still have, however, to establish a sign 
convention. We decide that 8 is to be taken positive if a 
turn about the axis of is from % towards This would 
make 

Yi - ft = -0(&-a*); y*-& = *(&-«i) < 4 ') 
with 4=4 = 0, lj = L Hence 

Yi — Pi = *fJCtn40 (Pm ~ ^m) (48), 

If we write 

l k 8 = m ki 

Yi - Pi = (Pm - O (49), 

For instance, if a m = 0, we have 

Yi - Pi » %A - *kP*l Yt - A * «0»A - »iA; 

y* - A = ®A - ®* A (60). 

For finite rotations we return to (44) or (41), keeping the 
positive sign in the second term. We may transfer the 
origin to a ( to save writing. The coefficient of JB* in y t — p i 
id 

** = (1 - cos 0) 44 - sin 6e ihm l m (51), 
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for k # i; if k = i, the coefficient is — (1 — cos 6) (1 — 44)* 
Thus 

“ (1 — cos 6} (44 — S«) — sin 0 c ik (52), 

where c ik is the antisymmetrical tensor corresponding to 
l m> namely 

/ 0 4 - 4\ 

(-4 o 4 (53). 

\ 4 -4 0/ 

Thus the displacement is represented in general by 
where b ik is a tensor of the second order, expressed as the 
sum of symmetrical and antisymmetrical tensors. The anti¬ 
symmetrical part is seen to be of the first order of magnitude 
in 8 and the symmetrical part of the second order* 

EXAMPLES 

1. Given that the general quadric surface is 

S = $A ik x t x k + r + C = 0, 

prove that the locus of the mid-points of parallel chords is a plane, 
and find the condition for tliis plane to be perpendicular to the 
chords. 

2. Find the condition that the line 

-tpc 0 


may touch the quadric S , 









CHAPTER III 


PARTICLE DYNAMICS 

The essence of particle dynamics is that the second de¬ 
rivatives of the coordinates of a particle with regard to the 
time are equal to functions of the position and velocity of 
the particle with reference to neighbouring particles. The 
relations therefore provide a set of differential equations 
to determine the coordinates. The equations for any 
particle can be put in the form 

m£i = EX, (1), 

where m is the mass of the particle, X, are the forces due 
to the other particles, and the summation is for aH the 
other particles*. It is a matter of experiment that this 
form holds when the axes are a certain type of Cartesian 
axes, which we call non-rotating and unaccelerated, or in 
brief dynamical 

If we take a different set of non-rotating axes at the 
same origin we have 

*f = (2)j 

and since the axes are not rotating the a i$ arc constants. 
Hence by differentiation 

® <Mi (3)j 

x/ - a ti x t (4), 

and therefor© velocity and acceleration are vectors. 

The force X, on the particle due to some other particle 
is measured by the contribution to mx 4 due to the other 
particle; that is, the part of mx { that would disappear if 
the other particle was removed to an infinite distance. It 

* Of. Jeffreys, Scientific Inference , Chapter vni, for a fuller analysis. 
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is actually measured by the acceleration. If we have the 
three acceleration components of the particle 1 due to the 
particle 2, we can find the acceleration in the direction of 
x / by the formula (4), and the form (!) will become 

mx/ = EX/ (5), 

provided that we define X/ by the rule 

X/ = a,jX, (6). 

The meaning of the force in any direction requires de¬ 
finition in any case: if we define it by its relation to the 
acceleration component in that direction, it follows auto¬ 
matically by (6) that force is a vector. It follows that the 
sum of any number of forces, obtained by adding the 
components separately, is also a vector. But the practical 
importance of the idea of force arises equally from the fact 
that in many cases the force components are known once 
for all from experience as functions of the coordinates and 
velocities. The form (6) then says, as a general principle, 
that the forces in any direction are additive. 

If we form the contracted or scalar product of ?hx, and 
EX f by the vector x iy we get 

mx i x t = EA^x, (7), 

d 

the left side of which is (£mx/), {The square of course 

implies the product of x, by x, and therefore the summation 
for i = 1,2, 3.) We write 

T = = \m (xj* + x 2 + x 3 3 ) (8), 

and call T the kinetic energy of the particle. Then by in¬ 
tegration with regard to the time from t 0 to ^ we get 

=eJ XfXfdt (9). 

But in any short interval of time dt , x t dt is the increase 
of x it namely dx t . Hence the right side is equal to 

3 JCT 
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21 XidXf taken from the initial to the final position of the 

particle. We call the scalar Xidx { the work done on the 
particle by the force X t in a small displacement dx £ . Then 

2 | X t dx f is the total work done by all the forces on the 

particle during the motion. We have therefore the scalar 
relation that in any motion of a particle 

Increase of kinetic energy = Work done on the particle 

( 10 ), 

In a system of particles we may add up the equations ( 10 ) 
for all the particles. We now take T , the kinetic energy of 
the whole system, as the sum of those of all the particles. 
Then we get 

[V] * = |s J* X t dxt (II), 

the double summation implying summation for all pairs of 
particles; the first summation is for the particles producing 
the forces X< and the second for the particles acted on. It 
may happen that the X t are all functions of the coordinates 
alone, and not of their velocities, and that provided the 
initial and final values of the x { are the same the integral is 
the same however the x* may vary in the interval. If so, 
the integral is the difference between the values of a certain 
function U for the initial and final positions of the system; 
we call the system conservative and U the work-function* 
Then (11) becomes 

Ml-K < 12 >' 

or T — U = constant (13), 

This is the equation of conservation of energy. The quanti¬ 
ties on both sides are scalars, — U is often denoted by V 
and called the potential energy. 
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If U exists, then in all possible small displacements of 
the system 

22 Xfdxi = dU (14), 

where the dx f , Zn in number for n particles, are all in¬ 
dependent. The coefficient of dx t for any one particle is 
2X,-, the total force in the direction of x £ on that particle. 
Hence the force on a particle in the direction of x ( is 
dUjdXi, which is of course a vector. 

We can also take the scalar product of ( 1 ) by any set of 
small quantities S x t whatever that constitute a vector. 
Then 


mXiSxt = 2X t Sx f (15) 

is a scalar relation; but as the Sx £ are arbitrary we can 
equate their coefficients and regenerate the equation ( 1 ). 
If we now integrate with regard to t from (q to we get 


I mxihxidt = f 2X,Sx<dl (16). 

Ji* Jt» 

The left side, on integration by parts, gives 
f l* 1 f (l d 

|^x<Sx,J — J mXi ^Sxidt (17), 

But we can consider the x ( 4 - Bx ( as coordinates of a par¬ 
ticle in a motion differing slightly from the actual one; that 
is, at a given value of t, the coordinates are x t + Sx ( instead 
of x<. Then 

d n, d . k . dr 

di Sx ‘=dt (x ‘ + SXi) -dt x ‘ 

= Xi + Sx* — x t 


(18), 

since x* + Sx t is simply the varied velocity or rate of change 
of the varied coordinate x £ -f 8x f * Then 


ffi d rtj 

J mXi ^ 8a? * dt = J ?nx i Sx, dt 

= \ h 8 (Imxf) dt + O (Sit)* (19). 
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Then for every particle, to the first order, 

J‘* {S {\rnxf) + £A r ,Sa;,} dt = [ma,S*,J* (20). 

We ean now add up these equations for all the particles of 
the system. If V exists we can express the result in the 
form 

j h S {T 4- 17) dt = 2 (21), 

or, since the limits are not varied, we can move the S out¬ 
side the integral. If the variation is such that the initial 
and final positions of the system are unaltered, Sar ( = 0 
when i = (u or , and we have 

S [‘‘(T + V) rff = 0 
Jfa 

to the first order in the variations of the path. This is 
Hamilton’s principle. 


EXAMPLES 


1. If the Cartesian coordinates of every particle of a system are 
known functions of a set of generalized coordinates q f9 prove that 


dx i 



Bq ft 


SX.&c^SX.-^S q„ 

where the summation convention is also understood on the right 
side. Deduce Lagrange's equations for the case where the q r are 
all independent. 


2, The equations of motion of a particle are ■= — k± i 

where k Is constant. Prove that 

2T=- x i X i + ~ (lm£*) + ± (\k£?). 

Hence show that for a system in periodic motion, or in one slowly 
changing its state, on an average over a long time, 2T ~ X 4 , 
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DYNAMICS OF RIGID BODIES 

A rigid body is one such that whenever it is displaced the 
distance between any two particles of it is unaltered. Since 
three particles A,B $ C are in a straight line if the sum of two 
of the distances AB r BO , CA is equal to the third, it follows 
that straight lines are unaltered by displacements of a rigid 
body. Since when A , B f C are not in a straight line the angles 
of the triangle ABG are determinate functions of the three 
sides, it follows that all angles are unaltered by displace¬ 
ments of a rigid body; If three lines meeting at a point and 
fixed in the body are mutually perpendicular before dis¬ 
placement, they are still perpendicular after the displace¬ 
ment. 

The equations of dynamics in the form mXi = X i are 
true with respect to dynamical (that is, non-rotating and 
unaccelerated) axes. Let any particle Q of a body have 
coordinates x t with reference to dynamical axes at 0, Then 
let the body be displaced in any way, and let the particle 
have the new coordinates x t \ We require the relation 
between x/ and x 0 Consider a particle of the body, P say, 
whose coordinates before and after the displacement are 
iit and a /. Put 

x t = a f + y { ; %/ = a/ + fft (1)* 

Then y ( and y/ are the coordinates of Q with respect to 
axes at P parallel to the dynamical axes before and after 
the displacement. Also if we imagine the original axes at 
P to be specified by the particles on them, these particles 
in the new position still specify a set of rectangular axes, 
with respect to which the coordinates of Q are still y*. If 
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the cosine of the angle between the y k axis in its new posi¬ 
tion and the y £ axis in its old one is a ik , we have therefore 

Vi « ^tkVt ( 2 )> 

whence x/ = a / + a ik (a:* — a k ) (3). 


The displacement of Q is 

a/ - a, + a ik y k - y t (4), 


We can prove that there is a straight line of particles in the 
body such that the displacement of Q is the same as that 
of A For such a particle we should have 

a ik y k - y t = (a ik - S rt ) y k = 0 (5), 


and these three homogeneous equations in y k are consistent 
provided 

I ! £% “ 8^ 1 1 = 0 (6), 


that is, 


a x i — * 


«!! 


,= 0 

<h* 


= 0 


(7). 


£% a _ 1 #23 

£*33 — 1 

Now the determinant || a ik || is unity, and each element in 
it is equal to its first minor. If we expand (7) we get 


lla^H — {(&u Q 22 “ ■+■ (£*22^33 — # 23 ^ 32 ) 

+ (U33 £?ix “ £*3(^13)} "b (®n "I" d* G33) “ 1 == ® (®)i 

since the terms all cancel. Hence (5) have an infinite 
number of solutions, all proportional, and the points there¬ 
fore lie on a straight line. 

If we take any two planes through this line, the angle 
between them is the same alter displacement as before, 
and therefore all planes are rotated through the same angle. 
Thus any displacement of a rigid body is equivalent to a 
displacement of a particle of it combined with a rotation 
about a line through that particle. If the angle of rotation 
is 9, we have, by comparing (5) with (52) of Chapter 11 , 


&ik ~ a ih "" 


8 * 


- {Ith - (1 - cos 0) - € ikm l m sin 9 (9), 
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where l t are the direction cosines of the axis of rotation. 
Thus the a ik are determined in terms of this axis and the 
angle of rotation. Conversely, if we take the symmetrical 
and antisymmetrical parts, 

(M* — 8 *fc) (I — cos 0 ) = | (b ik + 6 fcl ) ( 10 ), 

i fafciwi — €,*m) l m si * 1 ^ } (bik ~ b ki ) (H), 

But for a given m f with ikm in cyclic order, 

€ *im “ € ikm = ~ 2 (12), 

Thus l m sin 9 = - } (b £t - b ki ) (13), 

and 

sin- 9 = l {(612 — & 2 i) s + (^23 — ^a) a + (831 “ ^is) 2 } 04). 

Now suppose that the displacement is small. Then a ik 
have nearly their values for zero displacement, that is, 
S fJt , The direction cosines of the y k and y r axes in their new 
positions are a ik and a it , and thus, if k ^ l , 

a ik a n = 0 (15), 

and, if k = l, *tk a n = I ( 10 ), 

But in (15) for i = k t a ik is nearly 1 and a it small; for i = i, 
a n is nearly 1 and a ik small; for i not equal to k or l f both 
a ik and a u are small. Hence to the first order 

4- & k i = 0 (17), 

If k = i, then for i ^ k or l both terms of (16) are small of the 
second order, and therefore for i = k or Z, a ik = 1 4 - a 
second order quantity. Thus to the first order 

a ik = — b ik (18), 

where b ik is an antisymmetrical tensor. 

The displacement of Q is 8 4- a ik y k — y u where ha { is 
the displacement of P 7 and is a first order small quantity. 
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Now let X i be the total force acting on the particle at Q. 
In the displacement the work it does is the scalar 

BW = X ( (8a, + a ik y k - Vi ) 

= Xi {Sa, — b ik y k ) 

— XjScii — £5** (Xiy k + X k y f ) — lb ik (X { y k — X k y f ) 

(19). 

When i and k are interchanged, b ik and X t y k — X k y t are 
reversed in sign, while X 4 y k + X k y f is unaltered. Hence 

&W = XMi + | (y t X k - y k X 4 ) b ik (20), 
Now write y t X k - y k X t = L ik (21), 

so that L& is an antisymmetrical tensor. The second term 
in (20) is the sum of nine terms, of which three are always 
zero and the others equal in pairs. If wo replace b ik and 
L tk by the associated vectors, we have 

m = Xi5ai + L m b m ( 22 ), 

where the Sa, and b in are all independent of one another 
and the same for all particles of the body. Hence if we add 
for all particles of the body 

8 IF = (SA*) 8Ui 4- (2Z m ) b m (23), 

and the total work done in any given small displacement is 
determinate if we sum up the forces acting by the six 
expressions 2 X it Further, the contributions to 

2A%, from the internal reactions are zero. This follows 
at once if these reactions consist of equal and opposite 
forees between pairs of particles along the line joining 
them, and also has the justification that it leads to correct 
results. Then we may restrict to the contribu¬ 

tions from the external forces. This is d'Alemberts Prin¬ 
ciple. 

In the limiting case of continuous motion, we may con- 
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sidor the displacements that take place in a short time 
interval S£, and put in the limit 

Sa* = Ui$t; = a> ik 8t (24), 

and call u ( the velocity of P and oi ik the angular velocity 
of the body. Then the velocity of Q is 


= w ik y k 


(25). 


Now consider the centroid G f with coordinates de¬ 
fined by 

(2m) — EmZi (26)* 

It is usually assumed without proof that G is fixed in the 
body, though this is not obvious. But suppose that the 
particle at G with coordinates in the original position of 
the body goes to G r when the body is displaced, its new 
coordinates are given by 


{.r* — a t ) (27), 

and the coordinates of the new centroid G” are given by 
(2?n) W t 9 == Hmx/ (28). 

We have to show that G f and G” coincide* We have 


(2m) = 2m {a/ + a ik (x k - a k )} (2D), 

and therefore 


(2m) — #/) = Sro {a/ + a ik (x k — a*.)} 

- (2m) {V + && ($ k - %)} 

- 'Lma ik x k — Zma ik x k 
“ 0 (30). 

Thus the particle originally at G is displaced to the new 
position of the centroid, and therefore the centroid is fixed 
in the body. 

Now return to the equations of motion of a particle of 
the body 

mXi = X, 


(31). 
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By addition we form from these the equations 

Emi ( = EX, (32), 

and by cross-multiplication followed by addition for all 
the particles of the body, 

Em (x £ x k - x k x £ ) = £ fox* — o^X,) (3S). 

On the right sides of (32) and (33) the contributions from 
the internal reactions are zero, by d’Alembert 3 s principle* 
Also 

Emx t = Ema; # = ^ (Em) (34), 

where M = Em (35), 

the total mass of the body* Also if 

ha = Em (x t x k — $*£*) (36), 

we can reduce (32) and (33) to 

M$ t - EX, (37), 

| = 2£ tt (38). 

These are the fundamental equations of rigid dynamics. 
The three independent h ik are expressible in terms of the 
three independent co ik by (36), and we have therefore six 
differential equations for the x t and which determine 
them in terms of the initial conditions if the external forces 
are given. The motion of the body is therefore determinate. 
The principle of virtual work follows immediately. For 
if a body is initially at rest, x £ and h ik are zero, and the 
condition that they may remain zero is that EX, and EL ik 
shall vanish. But this implies, by (23), that in any small 
displacement of the body the work done by the external 
forces is a small quantity of the second order in the dis^ 
placements. Conversely, if there are six independent 
possible small displacements such that the external forces 
do no work, to the first order, in any of them, the co- 
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efficients of the independent 8 a £ and b m of (23) must all 
vanish* But these coefficients are the EX* and EL mf and 
therefore the x { and h ik do not vary with the time* The 
vanishing of the work done by the external forces in six 
independent small displacements is therefore a necessary 
and sufficient condition for equilibrium. 

The equations (38) may be put in another form. If we 
consider any moving origin 0\ not necessarily fixed in the 
body, with coordinates a u we can write for the coordinates 
of a particle with respect to 0\ 

y i ^x i - a t (39). 

If any vector associated with the particle, such as its 
velocity, momentum, or acceleration, or the force acting 
on it, has components u Sf we may form the an Asymmetrical 
tensor y t u k — y k u £ and call it the moment of the vector 
about O'* From (32) and (33) we can form the equations 

Em (a?,-#* — — (a £ Emx k — a k Emx £ ) 

- 2 [x 4 X k - x k X £ ) - (a £ ZX k - a k EX ( ) ( 40 ), 

where only the external forces make any contribution to 
the right side. But by (39) this is equivalent to 

'Lm ( y t x t - y k £ t ) = 2 {y,X k - y k X t ) (41). 

Therefore the sum of the moments of the mass-accolera- 
tion products about awy origin is equal to the sum of the 
moments of the external forces about that origin* 

If we denote the moment of momentum, or angular 
momentum, about 0 f by h ik , we have 

hx' = fait - y k £i) (42), 

and 

j 

-■ hit = 2i» (! ii£ k - y*Xi) + (i/i £ k ~ VkXi) 


( 43 ). 
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The second term may be written 

Em {(^ — it,) x k — (#* — a*) = — Em {a*#* — £*£*) 

4 =- M (iaA - d k x t ) (44), 

Thus (4i) is equivalent to 

^ ^ik + M ” &***) = E {yiX k — y k X t ) (45), 

In many important cases the second tensor on the left 
vanishes identically. This is clearly true if O' is fixed, when 
the time derivatives of the a, are zero; when the centroid 
is fixed; and, if both are moving, if the vectors and 
are parallel, that is, if the velocity of the moving origin 0 r 
is parallel to that of the centroid. The most important case 
is where the moving origin is identical with the centroid, 
when the last condition holds automatically. These terms 
also disappear if the moving origin is an instantaneous 
centre of rotation always at the same distance from the 
centroid. They vanish for a sphere or circular cylinder 
rolling down an inclined plane, but not for a rolling elliptic 
cylinder. 

If the moving origin is the centroid, 

A*' = {Vi (4 + y k ) - y k (*, + yi)} 

= iViVk - VkVi) ( 46 ), 

since 'Lmy i = Emy* — 0, 

by the definition of the centroid. Thus the angular mo¬ 
menta about the centroid are expressed completely in 
terms of positions and motions relative to the centroid, and 
the formulae for them have the same form as those for the 
total angular momenta with reference to a fixed origin, 

By (25) y t = - 0 J im y m (47), 

m being here a dummy suffix. Substituting in (46) we have 
hk = - Em K m y,y m - <*> im y k y m ) 

= B km W im - BtmWkm (43), 
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w r here B im is the symmetrica! tensor defined by 

B im = Em y t y m (49), 

and depends only on the masses and positions of the par* 
tides of the body. In the three-dimensional case, with i 
and k unequal, m must be equal in turn to i, k $ and the 
other value. Then, for instance, 

= (^21^11 H“ ^22^12 4“ -^23^13) 

' (^11^21 4 - 4 “ ^13^23) 

= {£>22 4“ B n ) -^13^23 (50). 

If we replace the antisymmetrical tensors by the associated 
vectors tills takes the foi'm 

^3 ™ ( ^11 + B22) <^3 -^ 13^1 

= (U U 4 - B 2Z -h £33) U)$ — -f- 4 “ ^33^3) ( 5 l)j 

and in general 

1=1 d (m o*^ (52), 

where A im is the symmetrical tensor given by 

— B im 

= (Emy* 2 ) S im 3jmy t y m (53). 

It evidently corresponds to the system of moments and 
products of inertia given in ordinar}' treatises on dynamics; 
A n , Am> are the ordinary moments of inertia A, B, £7, 
but A n , i4 3lJ A l2 are equal and opposite to the ordinary 
products of inertia F t 6 , //. 

It should be noticed that this reduction is characteristic 
of three dimensions; in a higher number of dimensions 
there is no analogous simplification of the form (48), 

The equations of motion of a rigid body then take the 
form 


MS, = 'LX, 

~h a ’ = XL a '; or^'-Z LJ 


(54) , 

(55) , 
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where L ik ' is the moment of the external forces about the 
centroid. 

It may happen that one point of the body is fixed. In 
that case we may take the origin at that point. Then in 
addition to the known external forces there is a reaction 
at the origin which can if required be found from (37); but 
the reaction has no moment about the origin and the 


motion is given by (38), But in this case 

Xi= — co im x m (5G)j 

and we find, by a process analogous to the last one, 

h ik = D kyn a> im — D im ti) km (57), 

where D fm ^ Ema '. t x m (58), 

In terms of the associated vector, 

h m “ O im W t (59), 

where G im - S im - (60), 


The correspond to the moments and products of inertia 

about the origin. They can be expressed in terms of those 
about the centroid; for 

~ (3?* + S tTn —* Erra + yi) (3^ + y r/i ) 

= (Em** 2 ) S im - ’SmW i x m 
+ A im (61), 

the terms linear in the y *s vanishing by the definition of the 
centroid. 

The relevant equations of motion then take the form 
J,K~-ZL m (62), 


where the L m are the moments of the externa! forces about 
the origin, the reaction at the origin making no contri¬ 
bution. 
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The kinetic energy of the body is given by 
2 T = Emi\ s 

= Em - rn im y m )* 

— 2mi ( a + Hmy k y m a> ik ti> m 
= Mxf + B km (63). 


In three dimensions w ik u> im is ze to unless i is different from 
both k and m. If i = 1, k = 2, m = 3, or if i = I, k = 3, 
m = 2 } Cf>u4) is = — o>,«j 3 . If i = I, k = m == 2, w n z = 
and has coefficient Thus w, 2 enters with & coefficient 
B 22 + B n , or A x . But w 8 to 3 has a coefficient — 2 B.& or 
A a + A 3t . Thus in all 

2 T = 3f.t* + A km w k u) m (64). 

When a point of the body is kept fixed, 

2T = C km w k oi m (65). 

We notice that the linear and angular momenta take the 
forms 

dT 3 T 

Mx 4 = w . ; h t * — (x t constant), 


h t — 


dr 

dw t 


m 


for origin at a point of the body held fixed. 

In one respect the forms (55) and (02) arc inconvenient. 
They involve the tensors A im and G im> which depend on 
the x t and and therefore in general change as the body 
rotates. It is more convenient to use such axes that the 
relevant tensor in the particular problem is constant. To 
achieve this the axes must rotate, and then are no longer 
dynamical axes. Suppose then that we have a set of axes 
x/ rotating in any way, and that their direction cosines 
with respect to the dynamical axes a?* are a u . All the usual 
tensor relations hold for transformations from the x/ 
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system to the x t system and vice versa; but the a iiy while 
specified for each instant, are now functions of the time. 
If Ui is any vector, we have 


u t = a^Uii %i = a u u / 
Then is another vector; and 


dUi du / 

1U ~ * 


_ 4 _ w 


(67). 


( 68 ). 


The component of this vector in the x/ direction is 

duj f da £4 , A , 

a u a a —fa + Gnuj { 69 }* 


The first term of this is 


5 du/ du/ 
bli ~dt - ~df 


(70). 


da 

Also is the x t velocity of a point at unit distance along 

the x/ axis, that is, a point with coordinates a a (i — 1, 2, 3) 
with reference to the dynamical axes. This velocity is 
— & ik to kU where 8 ik is the antisymmetrical tensor expressing 
the rotation of the rigid frame consisting of the moving 
axes* Then 


t daj * t A 

toiiUf ^ - — — a it Uf 8 (k a ki 

= * 8 l /u/ 

and the required component is 
du/ 


dt 


— 8 t /u/ 


(71), 


(72). 


If we use instead of 8 S / the associated vector, the three 
c om ponents b ecom e 

Oh — 8$ + @ 2 , u/ V®/ 4- u/8 / s 

< - u/8/ + u/0/) (73)* 
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If for instance the components of a displacement along the 
x/ axes are u/ t the formula (73) gives the components of 
the velocity along these axes; if u/ are the components of 
velocity, (73) gives the components of acceleration; if u/ 
are the components of angular momentum, (73) gives the 
components of their rates of change with reference to 
dynamical axes, and these are equal to the components of 
the moments of the forces acting. 

EXAMPLE 

Prove that 

- 


1CT 


4 










CHAPTER V 


EQUIVALENCE OF SYSTEMS 
OF FORCES 

An external force X £ acting at a point x t of a rigid body 
produces dynamical effects summed up in the vector Xt 
and the antisymmetrical tensor x t X k — x k X it A force lias 
a line of action; that is, if we take its resultant R given by 

R* = Xt* (I) 

we can define a direction Z< by 

X t - Rl t (2). 

By convention R is always taken positive. If a force X, 
acts at the point x t 4 rl u where r is arbitrary, we have 

(Xf 4- rZf) At Hh rZt) X £ = x £ X k x^Xf (3), 

Thus the dynamical effects are the same if the force X f 
acts at any point of a line through x t with direction cosines 
proportional to X,, This line is called the line of action of 
the force, and the force can be said to act along it. 


If a force has magnitude R and acts at x £ in the direction 
Z,, we have 


) 

Xi = Sl t 

(4). 


Ltk = R { x i^k 

(5), 

or 

% 

J 

BJ 

II 

(6), 


in the notation of the coordinates of a line. Thus the X, 
and L m are the products of the resultant into the six 
coordinates of the line of action. 

The moment of the force X, about a point is 

(#i — a t ) X k — (x k — a*) X, = R (l ik — a £ l k + (7). 
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If we take the associated vector, its component in the 
direction m m is 

R (m m lj ~ = R 4- Umf) (S), 

where m/ are the other coordinates of the line through a t 
in the direction m,. We notice that 

m 4 l / 4- hm/ = d sin 9 (9), 

where 9 is the angle between the two lines and d the length 
of their common perpendicular. We may call Rd sin 9 the 
moment of the force about the line (m u m/). 

By d'Alembert's principle, the motion of a body is un¬ 
altered if to the forces acting on it we add two equal and 
opposite forces acting at the same point, or, by (3), along 
the same line of action. 

Now consider a pair of equal and opposite forces X f and 
— X f acting at points a t and b € * They clearly make no 
contribution to EX*. Their contribution to HL ik is 

(a £ X k — a k X ( ) — (b i X k — ftjtXf) 

= (a< - b 4 ) X k - (a k - b k ) Xf (10). 

Since a t and b £ are equally affected by any motion of the 
origin, the contribution of such a pair of forces applied to 
definite particles to both EX* and Ei lJt is independent of 
the position of the origin. Such a pair is called a couple , 
and its contribution to S L ik is called the moment of the 
couple . 

If the vectors a t — bi and X* are both perpendicular to 
a line with direction cosines n it the components of E L m are 
the components of a vector along this line; while the forces 
act in the same plane perpendicular to the line. This plane 
is called the plane of the couple, and the line an axis of the 
couple. Evidently equal couples in parallel planes are 
equivalent. The magnitude of this vector is Rd, where 
R 2 = Xf- and d is the perpendicular distance of a { from a 
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line through 6, parallel to X ( ; and its components can 
accordingly be written Rdn m . 

Any system of forces is equivalent to a force at an 
arbitrary point together with a couple. For with each 
Xf acting at x { we can associate a pair of forces ± X, at a t . 
Then the system is equivalent to 2X, at together with 
a set of couples whose total moment is 

M ik - 2 {(*< - a t ) X k - (x k - a k ) X t ) (11). 

But this is equivalent to a single couple; for we have only 
to make n m proportional to the M m and Rd equal to their 
resultant. 

If L ik are the moments about the origin, 

Afr = 2 ~ x k^i) (I2)j 

and if for brevity we replace 2X* by simply X i$ 

^ik = A* ~ (&iXt — (13), 

or M m = L m - € ikm a { X k (14). 

Evidently X { 2 is a scalar and independent of a { * Also 

X m M m = -An Am “ € iJem a tX k X m (I®)- 

The first term is a scalar and independent of a t . The second 
is identically zero. Hence X t 2 and XfMi are scalar in¬ 
variants. 

The system is equivalent to a force X ( and a couple M t 
at These vectors are parallel if 

An — *ikm a iX k — pX m ( 16 )j 

where p is a scalar length called the pitch. These give three 
linear relations between the three a { and the pitch, and we 
therefore expect a single infinity of solutions. But if we 
take the scalar product of (16) by X m we have 

X m L m = pX m * 


(17), 
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so that p is the same for all admissible values of With 
this value of p, (16) represents three planes with a line in 
common. If we change a t * to a* + &X i? where a is a scalar, 
the left side of (16) is increased by 

tmtoXtXk = 0 ( 18 ), 

and therefore if a t is one point on the line, all points on the 
line through a * parallel to Xi satisfy the conditions. This 
line is the central axis of the system; the system is equiva¬ 
lent to a force R along the central axis and a couple 0 
about it. This expresses the system as a wrench , If we take 
the central axis as one of the coordinate axes, we have, 
since X t 2 and X t M { are scalars, 

R 2 - X t 2 (19), 

GR+XtMt-pXf (20), 

and therefore G , 2?, and p arc determined. 

The system can also be reduced to a couple parallel to 
a preassigned plane together with a force. For if S is the 
couple, and % are the direction cosines of the normal to the 
plane, and if the force acts through a it we have, for the 
moments about a it 

M m = L m — c fkm &{ X k (-1), 

and also M m — Sn m (22). 

We have three equations to determine the a t and S, Again 
there are a single infinity of solutions. But if we take the 
scalar product by X m , we get 

A m = X m Gikm&t At X m 

- OR (23), 

so that S is determined provided n m X m is not zero, that is, 
provided the resultant force is not parallel to the plane. 
Then the equations 


An — ^ikm ■Als — Sll m 


(24) 
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determine a line parallel to X u which is the line of action 
of the force. 

If l t are the direction cosines of a line through a it the 
moment of the system about this line is 

“ A An 

= l m L m + 4'X* (26), 

where Iff are the other coordinates of the line. If this 
moment vanishes the line is called a null line of the system. 
If 6f is another point on it, 

An (Aj “ a m) = €ikm a iX k (b m — & m ) (27), 

which shows that b { Ues in a definite plane through a i . All 
null lines through a point therefore he in one plane. This 
plane is called the null plane of the point. 

All null lines in a plane pass through a point. For the 
system can in general be reduced to a couple in the plane 
and a force whose line of action intersects the plane in one 
point. Then the system has no moment about any line in 
the plane through this point, which is the null point of the 
plane. 

Any system is equivalent to two forces, one of which can 
be made to act along a given line. For let the lines of 
action pass through and b ( , and have direction cosines 
and m ir and let the magnitudes of the forces be 8 and T. 
Then we have six equations, 

X* - 8li + Tm f (28), 

L. m - € {ktn (atSljt + b t Tm k ) 

- SIJ + TmJ (29). 

The coordinates of the first line being given, these arc six 
linear equations to determine the six coordinates of the 
second line and jS and T r But we have also 

nit 2 = l; — 0 


(30), 
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so that we have in general just enough equations. Two such 
lines are conjugate lines . Clearly any hue intersecting two 
conjugate lines is a null line. 

It can bo shown easily that 

S {XJ/ -f LM ~ L i X i (31), 

r* = X, a - 2£4X< + s* (32), 

whence (28) and (29) determine the coordinates of the 
second hn© explicitly. 


EXAMPLES 

1. A system of forces is reduced to a force at P together with 
a couple; P is chosen so that the couple is parallel to a given 
plane. Show that the locus of P is a straight line parallel to the 
central axis. 


2. Show that 


X? 

% 


is a point on the central axis, 

3. Two systems of forces are given by (X i9 L i ),(Y i9 X t ). Show 
that X i K i + is invariant. 
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CONTINUOUS SYSTEMS 


In problems involving volume and surface integrals we find 
it convenient to denote elements of volume and surface by 
dr and dS respectively; both are always taken positive. 
They are of course scalars, The direction cosines of the 
normal to an element of surface are usually denoted by £*; 
in most cases the normal is drawn otriwards from the region 
under consideration. 

Green’s Lemma tabes the form 


where V is a scalar* 


where 


V* = 


f Si *-JJ 

( 1 ). 

diary, if 


U f — dVfdx t 

(2), 


(8), 

3* 3« 3 a 9* 

(4), 

9* ( a — 3*1 2 3*2 s 3*3 3 


and 3 jdn denotes differentiation along the outward normal. 
Stokes’s Theorem takes the form 


J.-MJ -£)+<■(£-£) 

- | J h**m duJdXk dS 




(5). 


The integral on the left is round a closed contour C , On 
the right is the direction cosine of the norma! to any 
element dS of a surface S that fills up the contour; the 
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integral is over S; and the sense of the normal is such that 
if the contour is described in the positive sense about any 
axis * l{ is taken positive when the normal is in the direc¬ 
tion of x t increasing. 

The gravitational potential of a distribution of matter is 
given by 

V = ( 6 ), 

where m is an element of mass, r is its distance from the 
point where F is to be found* and / is a constant equal to 
6*66 x 10" s when m and r are measured in grams and centi¬ 
metres and the time in seconds. When the mass is distributed 
continuously over a surface or through a volume, m must be 
replaced by adS or pdr respectively, where tr and p are called 
the surface and volume densities. The work in displacing 
a mass m* through a small distance is m'd F, 

The electrostatic potential of a system of point charges 
is given by 

V = fZ e r ( 7 ). 

where e is a typical element of charge. If e is in electro¬ 
static units of charge and r in centimetres*/is + I, The work 
in a small displacement of a charge e' in the field is — e*dV\ 
The usual relations follow, that in free space in both 
eases 

V*F - 0 (8)* 

and in space occupied by matter of finite density 

V*F = — 47 rfp (9), 

Also we have Gauss’s Theorem 

JJ |£ dS = - 4rr/S'«t or - 4rr/S'e (10), 

where the summation is over all the masses or charges 
within the closed surface S . In crossing a surface where 
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there is a finite surface density a, d Vjdn has a finite dis¬ 
continuity — 4 nfa. 

All these relations are scalar in form. For a gravitational 
field the force on a small particle of mass tn' is 

m'X t = m'dVjdXi (11), 

while for an electrostatic field the force on a small charge 

e , X t ~-e , dV/dx i ( 12 ). 

The vectors are called the intensities of gravitational 
and electric force respectively. The analogy in form be¬ 
tween (6) and (7) is constructed for mathematical con¬ 
venience; the difference in sign in (11) and (12) embodies 
the physical difference that whereas two positive masses 
attract, two positive charges repel. 

In a gravitating system we may construct a work- 
function 

W = ( 13 ), 

where the summation is over all pairs of particles m pi m Qi 
and Tjx ia the distance between m p and The case where 
p = q is excluded. Then the force on the pth particle is 
3 lV/(dXi) p . The potential at m p due to the other particles is 
V v = S 'fwiJrpQt the accent indicating that the caae where 
q = p is excluded from the summation ; and the function 
V P = 2 W, since each pair of particles is counted twice 
in this double summation. The function l'Zm p V P therefore 
plays the part of a work-function. Similarly in electro¬ 
statics the function l^e P V p plays the part of the work- 
function with its sign reversed, that is, of a potential energy. 
These results may be generalized to the case of continuous 
distributions; thus we can replace these functions by 

W-llljpVdT+ljjaVdS 


( 14 ). 
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where the first integral is tlirough all regions where p is 
finite and the second over all surfaces where a is finite; the 
first integral may therefore be taken tlirough all space, 
excluding the surfaces where there is a surface density. It 
can be shown that with proper precautions about the de¬ 
finition of V the restriction p ^ q gives no trouble provided 
V is everywhere finite. 

Now if we consider the integral through all space except 
thin laminae surrounding the surfaces where a is finite* 



where the dn in the first integral is out from the region of 
integration and therefore towards the surface where a is 
finite. On the two sides of such a surface the values of V 
differ by an indefinitely small amount, and for the two 
dV 

sides together ^ dS = infadS, by Gauss’s Theorem. 

Hence the first integral in (15) is equal to — ^ jj gV dS 

taken over the surface and cancels the second term in (14). 
Hence 



through all space. In consequence of this form we may say 
that the gravitational work-function, or the electrostatic 
potential energy, is per unit volume, where R is the 

resultant of the appropriate intensity vector. 

When the properties of the medium vary from place to 
place, V is no longer of the form 2/e/r, and V 2 V is no longer 
zero. But a potential still exists; if a small charge e* is 
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moved from a point Ptoa point Q, the work done is still 
the same whatever the route, and may be denoted by 
e' (F P - V Q ). The treatment is suggested by the fact that 
two similar condensers with the plates at the same po¬ 
tentials, but with air between the plates of one and another 
material between those of the other, have charges in a ratio 
K depending only on the media. dVjdn on the outside of 
the condenser being small compared with its value be¬ 
tween the plates, we infer that the charge per unit area for 
the same distribution of V is related to the discontinuity in 
KdVjSn, where K depends on the material. This suggests 
in turn that Gauss’s Theorem must be replaced by 

jjz|?dS = -47r/2e (17), 

where the summation is for all charges inside S. Then 
applying this to the two sides of any surface we have 

[*s]—^ (18 >' 

and applying it to a region with a finite volume density we 
have 

£(*§£)*• ( 19 ), 

and as this must hold for all such regions. 



These equations are all homogeneous in V, />, c j; licncc 
the potential due to any set of charges is proportional to 
the charges if all are altered in the same ratio. Using this 
principle we can show by the usual method that the energy 
of a distribution is 

Ve — l \ |\pVdr + \ <jV dS (21)* 
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The first integral, applied to all space except thin regions to 
cut off surfaces where there are surface densities, gives 


-s yiKK-D"’ 

- - i* + £ JJJ*©'* < 22 >- 

The surface integral cancels the integral £ jjtrFdS; and 
hence 


W = 



m. 


so that the energy can be considered equal to KR 2 /Snf per 
unit volume* 

Magnetism may be treated similarly, starting from the 
assumption of volume and surface distributions of mag¬ 
netic pole strength, subject to the condition that the total 
pole strength in any solid is zero; or we may regard the 
ultimate magnetic unit as the doublet, which explains the 
need for the restriction involved in the former method of 
treatment. The potential at x t due to a doublet of strength 
M at the origin with its axis in the direction A* is 


V — yMXiZjr* (24), 

and if the doublet strength per unit volume in a solid is I 
in the direction A* we can introduce the intensity of mag¬ 
netization at $ it the vector IX it and say that the 
potential at x { is 

V = y\\\A i ^J i dr (25), 


where A i corresponds to the point £,• and dr = d^d^d, f s; 
y is a constant. The magnetic force in free space is 


a< = - 


cV 

8x t 


(26). 
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The potential is therefore equivalent to that due to a dis¬ 
tribution of magnetic poles l t A t per unit area over the 
boundary and — dA^dii per unit volume through the in¬ 
terior. 

Within a solid special treatment is needed. To define V 
or the force at x t , when x £ is within a solid, we must imagine 
a small cavity made about x f , the intensity of magnetiza¬ 
tion everywhere remaining as before* and consider V and 
a* within it; then the values of V and a ( at are defined to 
be the limits of those in the cavity when the dimensions of 
the cavity become indefinitely small. This process leads to 
little difficulty in gravitational and electrostatic problems, 
but in magnetism the limit of the force is found to depend 
on the shape and orientation of the cavity. The force in 
the cavity can be written 

*<--£ (28 >' 


where V is given by (27); in the first integral the normal is 
inwards towards the cavity. The contributions to X { from 
the volume integral and the outer boundary are of the same 
form as for gravitation, and give no trouble. If the cavity 
is a cylinder in the direction of the intensity of magnetiza¬ 
tion, IjA < is zero over the sides and equal to /, the resultant 
intensity of magnetization, on the ends. Such a surface 
density over the ends in the limit contributes nothing to V ; 
if the radius of the cylinder is small compared to its length 
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it also contributes nothing to X t ; but if the cylinder is 
of disc-like form it contributes 47 ry\ t l to X t , Hence if we 
take V for the complete body and define by 


a,- 


dV 

dXi 


( 20 ), 


a t is the value taken by X ( in a thin cylindrical cavity 
parallel to the intensity of magnetization. The force in a 
flat cylindrical cavity with its generators in this direction is 

a t = cii + 47ryA f / = a t + 4vyAi (30), 

Evidently and are both vectors; the former is called 
the magnetic force and the latter the magnetic induction. 
The theory of susceptibility and permeability may then be 
developed as usual* Also (27) shows that V is continuous 
across a boundary; but l f a t = BV/dn has a discontinuity 
— ^iryliAi ; whence Z*a* is continuous across a boundary. 
The mutual potential energy of two doublets M and M r 
at x t and x t \ oriented in directions A* and A/, is 

W = X ( 'M' ~ (31), 

where V is the potential at due to the magnet at x t ; this 
gives 

W-yMM-X! 

= yMM' {^‘ - 3 X * ~ ** Xi ' ~ * j (32) 


yMM' 

r* 


(cos e — 3 cos 0 cos O') 


( 33 ), 


where e is the angle between the axes of the magnets and 
6 and O' are the angles made by the axes with the line 
joining the centres. 
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If the second magnet is turned through a small angle Sip 
about a line with direction cosines n t , 

8 V = e ikm n k A m ' Sift (34), 

by (49) of Chapter n. Hence, by (32), 

sir 

sir = ~ SV 


= yMM' - 3 cos 6 


so that the couple about a line parallel to the axis of a:* is 
M k = - 3 cos (35). 

Hydrostatics and Classical Hydrodyiftamics * The internal 
reaction in a fluid across an element of surface dS is a 
pressure pdS normal to that surface* If the density is p, 
the bodily force per unit mass X ir and the velocity of the 
fluid at is u it the acceleration of the fluid is found, by 
considering a small parallelepiped, to be given by 



lx] + pXi 


(I). 


If u ( is given in the Eulerian way as a function of the 
coordinates x ( and the time £, the operator djdt , giving the 
rate of change of any element associated with a given 
particle of the fluid, is equivalent to 


d d t 3 

dTef + ti *er fc 


( 2 ). 


If we consider the circulation Q around any closed 
circuit C in the fluid and moving with it, defined by 


fi =* 


L 


UfdXf 


(3)i 
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we have 

= - f - fat + f X t dXi + [ u ( du t (4). 
JcpCZi Jc Jc 

The last integral is which always vanishes because 

when we move round the contour we come back to the 
same point, where the velocity has its original value. Also 
if X t is the gradient of a single-valued potential, as when 
the bodily forces are due to gravity (the commonest case), 

I X,dXi is the change of this potential round the contour and 

is zero. Again, if p is a function of p only, as in an incom¬ 
pressible liquid or a gas at uniform temperature, the first 
integral vanishes and 

S-o w. 

dt 

If then Q is ever zero around a circuit it remains so per¬ 
manently. This is true if the fluid is initially at rest and is 
set in motion by solids moving in it, and in various other 
cases of importance. But the vanishing of Q for all circuits 
is the condition for the existence of a velocity potential <f> 
such that 

u <=li (6) - 

In this case we can rewrite the equations of motion in 
the form 


cHtf 1 y 

si +u ‘m,--ps, +x ‘ 




and multiplying by dx ( and adding we have 


<8X 

3 Ui _ 3 2 <f> _ 3 % 


since 


(a). 
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— [ — -f- U + constant 

(10), 

P dU 

Ai ~d^ 

(ii); 


This shows that for all contemporaneous variations 
d<f> 
dt 

where 

u k 2 is the square of the resultant velocity q. The constant of 
integration is not necessarily the same at all instants and 
therefore may be a function of the time. Hence we have 
the Bernoulli integral 

d l + iq * = -j&+U + F(t) 

The rate of change of mass within a given small paral¬ 
lelepiped dx x dx 2 dx s is equal and opposite to the rate of 
outflow; hence gp t g 


( 12 ). 


-rdr= — ^ {pu^dr 
at ox £ r 

and we have the equation of continuity 
dp 9 , v 
dt ~ Sx t ^ 

dp 9 Uj 

OT dt = ~ p dx i 


(13), 

(H). 

(15). 


Vectors with given Divergence and CurL We sometimes 
have to find a vector u ( such that 


OX ^ 


( 1 ), 

( 2 ), 


du k diti _ 

0 x t dx k 

where A is a given scalar and o> ik a given antisymmetrical 
tensor. We want particular integrals of these equations. 


Evidently if 
where is any scalar, 


ds 

Ui = ~~ 
dXi 

V 2 <£ = A 


<3)> 

(4), 
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$ = - 

i\ j 


(G), 

where 

is to be evaluated 

at x t 

and A at 



dr = 


("), 

and r is 

the distance from 

ft to 

*c 


If 

dF m 

_ 3^ 

8F m 

(S). 


3a: m 

~* ikm dx k 

where F 

, is a vector such i 

that 




we have 
Then also 


=0 
ox t 

dtif _ d 2 F 

cx { ~ €ikm dx i Sx k 


” = 0 


(»)- 

( 10 ). 


du 


du 4 


“ ~dx t dxM “ fiJsn a* 


Su k 

a 


dF, 

*“ dx] ekps dx, 

3W. 

dx t dx v 

■ (3f 3 Smp— SfpSpu) 

3 *F. 


= & 


- 8 


&F t 
cXi dx p 

0*JP™ 


i8 dx i 9x 

a_ _ &K 

dx m BXi dx/ 


ip dx t dx p 



= - v»j- m 

(11) 

by (9). 

Thus if 

we shall have 

j 1 = _r f 

" 4l7, 

11T* 

(12), 


3w fc 3 k< , eht ( _ e 

dx t 3x k a}ik ’ Sx t 

(13), 
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provided (9) is satisfied. But 

bf. 




(14), 


since <u m is a function of $ { alone and r a function of 
x t — ii . Applying Green’s Theorem to all space except a 
small sphere about x { , we get 

dFr 
dx. 


= - j- Lim ff l ^dS + ~ Lim 
dx m 4 it JJ r 4 it 


^dr 
r o$ m 


(15), 


since F m has the form of a gravitation potential and 
cFnjdXf that of a gravitational force. But the first integral 
vanishes in the limit when the sphere becomes very small, 
and the integrand in the second is zero provided the com¬ 
ponents u ( exist, for 

du> m 3 (du k 0«A a (161 

W^d^Xdx, dxj ~dx,-dx m u { ’■ 

If then we are given a scalar A and a vector a> m such that 
its divergence is zero, a solution of (1) and (2) is 


w _ H,*?*_dj\ 

* dx f + dx k cx m 


(17), 


where </> and F m are given by (6) and (12) and ihn are 
in cyclic order. 

This analysis has two practical applications. In hydro¬ 
dynamics u t is the velocity and aj m is twice the vorticity, 
denoted by 2tj m in Chapter lx. Here the divergence of the 
velocity and the vorticity may be given through all space, 
and the velocity (17) satisfies the conditions. If the actual 
velocity is v t , we may put 

V t = Uf + u/ 

du/ = 0 . du/ du/ _ 

ox, ’ dx t vx k 


and then 


0; 


0 


(18), 

(19). 
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Thus u/ is the gradient of a scalar <f>' satisfying Laplace’s 
equation. There is no such scalar that makes the velocity 
finite everywhere, including at an infinite distance, except 
a linear function of the coordinates. Hence u/ is the same 
everywhere. 

If there are solid boundaries or free surfaces at a finite 
distance u/ may not be constant. If the region where 
vorticity is present does not extend to a boundary, there is 
no contribution to (14) from points outside this region; we 
therefore take (14) through a region large enough to contain 
the whole of the vorticity. Then the surface integral in (15) 
must also be taken over the boundary of this region, but 
still vanishes, and (17) is still a solution. But (17) may not 
satisfy the boundary conditions, and then we must add to 
U( an irrotational solution chosen to make the whole velocity 
satisfy them. 

In electromagnetism u ( may be the magnetic intensity 
and u> m the electric current across unit surface in a plane 
of x m constant. 

In many cases u, has no curl outside a limited region of 
very small cross-section. This is often true in the motion of 
a real fluid, when the region may be called a vortex fila¬ 
ment, and in magnetism, when the region is a wire carrying 
an electric current. The former statement may be expressed 
also by saying that the motion is irrotational outside the 
vortex filament; the latter says that magnetic forces due to 
electric currents have a potential. In either case the integral 

Q = f u ( dXf (20), 

Jo 

taken around a closed circuit, is zero if the circuit can be 
filled up by a surface not cutting the filament or the wire, 
and has a constant value for all circuits that cannot be so 
filled up. The two conditions are mutually exclusive, and 
therefore the critical region itself must be a closed circuit. 
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In each case A vanishes. Then <£ = 0 , and the component 
F m is given by 

'•-sJJJt* <21) 

taken through the critical region. But if we consider the 
contribution from an element between two planes sepa¬ 
rated by d$ m , and call the element of surface in a plane 
parallel to these dS, we have 

dr = dg m dS ( 22 ), 


[ = Cl 


(23), 


where the line integral is taken around the boundary of 
the filament. Hence 

f = R (iis 

m 4rrJ r 

the integral being taken around the length of the filament. 
Also 


(24), 


ra » («i 

1 4 t 7 (dx k Jr dx m J r } 

1 j Bx t (r) d ^ m 


£1 
: 4tt 


■£j“»^r )lmda (25} > 


where ds is an element of length of the filament and l m a 
direction cosine* 

In hydrodynamics Q is the circulation around the fila¬ 
ment* In electromagnetism the unit current is such that 
if it flows in a circle of radius 1 cm* it produces magnetic 
force at the centre* If we take the circle to be in the plane 
of constant, with its centre at the origin, we have 

ds = d $; r = 1 ; = z* x z = 0 ; 

= cos 0 , f t = sin 0 , — 0 ; l x = — sin 0 , 4 «= cos 0 , 4=0 
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and 

J *ikm = 0 for * = 1 or 2 

and *= 2 tt for t = 3. 

Also for unit current %, the magnetic force, = 2tt. 

Hence, in this case, by (25), 

Q ™ 477, 

and in general, if the current is /, Q = 4t tI s and 

I d8 (26) 

= (27). 

This may be transformed by Stokes’s Theorem into an 
integral over a surface with the wire as its boundary; thus 


L- ** 


dS 

d n - 


r f T 7 _ ® £k — 

u i = 1 jj gj" “^3 

■“ - 1 1[^» {S.mS.. — 8(.8 m ) (-) 

“ - 7 1 / L { S ‘-V (r) _ sSlZ (?)} ® 

=- 7 KU(-> 


dS 


(28). 


a 2 / 1 \ 

But — l m 57 —^— ( - ) is the magnetic force at x { due to 
cx t \rj 

a doublet of unit strength at with its axis in the direction 
l m * Hence the force at any point is equivalent to that due 
to a distribution of doublets over the closing surface, with 
intensity / per unit area and directed normally to the sur¬ 
face, Such a distribution constitutes a magnetic shell of 
strength L 
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It is customary to assert as a fundamental postulate that 
the magnetic force due to an electric current is equivalent 
to that of a magnetic shell, but I think this course un¬ 
desirable. The magnetic shell does not exist in nature and 
direct experimental test is therefore impossible. Further, 
though it gives the same force, it does not give the same 
potential; the magnetic potential due to a magnetic shell 
is a single-valued function with a finite discontinuity at the 
shell, while that due to an electric current is a cyclic 
function with no discontinuity except at the wire itself. 

Mutual Energy of Electric Circuit and Magnetic Field. 
A magnetic pole tn at x t is under a force mu ( ; but 



and — l m (-) is the magnetic force normal to the 

^ rr 3 /1\ 

surface at due to the pole. Hence — J | d& 


is the total flux across the closing surface, or through the 
circuit, of the magnetic force due to the pole. Also if V is 
the magnetic potential due to the current. 


Ui = - 


dv 

ox i 


<30), 


and therefore the correct forces are given by taking 


F = (31), 

where A T is the normal magnetic force due to a unit mag¬ 
netic pole. This can now be generalized; we say that in 
general the mutual potential W of a current and a set of 
magnetic poles is 


W = ~1 


NdS 


(32), 
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where N is now the normal magnetic force due to the whole 
of the rest of the field. This form may be further extended 
to express the mutual energy of two electric currents. 
From this equation we can infer, as is done in the standard 

works of Jeans and Livens, that if j J NdS varies with the 

time, the variation generates an induced e.m.f. — ^ [) NdS 
in the circuit. 

If we have two circuits carrying currents I and J, their 
mutual influence is expressed by the statement that their 
mutual energy is of the form (32), where N is taken to be 
the force near £< due to the current J ; points on the latter 
circuit can be taken to be given by x t . Then 

IF = IJM (33), 


where 


M 




(34), 


and u i is the magnetic force at f ( - due to unit current in the 
second circuit. But if is the vector potential due to such 
a current, 

dF m 


% — ; 


3fi 


M 


— j | 


,dS 




= J F m <Ur, 

taken around the circuit. But 


= J 


dx» 

7 


taken around the second circuit; and therefore 

‘dx m df m 


M 


=!! ! 


(35), 


(30) 


(37), 


(38), 


(30), 


taken around both circuits. This gives the required form 
for the coefficient of mutual induction of two circuits. 














CHAPTER VII 


ISOTROPIC TENSORS 

A tensor is called isotropic if its components retain the 
same values however the axes are rotated. We have 
already obtained three examples, namely S ( * f e iim , and 

There are no Isotropic tensors of the first order* For if 
u 4 was such a tensor, let us give the axes a small rotation 
expressed by the antisyinmetrical tensor c ik . Then in the 
new system 

tt/ — (Sjrj — Ctf) % = CtfUj (1), 

and this can be equal to u ( only if 

CtfUj = 0 ( 2 ) 

for all admissible values of the c iS . Thus 
%% 4- u 2 + %2/a = 0| 

CZL*H + C 22% + = 0? (3)* 

c 31 % 4 - c^u 2 4 = OJ 

But Cjj = — Cz 3 = 0 , while c 12l ,% arc independent and 

equal and opposite to the components obtained by inter¬ 
changing suffixes* Hence (3) can be satisfied only if 
= = ^ 

and therefore there is no isotropic tensor of the first order 
other than zero* 

11 u ik is an isotropic tensor of the second order, 
ik = (8* — (Sjtj — c kl ) u it 
— u ik — 

= ^ 1 *(^) 
to the first order, for all values of i and k, Hence 

°ii u ik + &ki U ii = 0 


( 2 ). 
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If i and k are unequal, take i = I, k = 2 * Since c n = Cm= 0 
we have 

^2^22 + C ia ft;i2 + %%1 H" ^23%3 =* 0 (3), 

and therefore 

^ ^ = 0; u 11 = u 22 (4). 

By symmetry u ik is therefore 0 if i ^ k, while u u = u n — 1 %* 
If i and k are both 1 , we have 

4 c 13 i% 4 612^1® + = 0 (5), 

wliieh is satisfied since every term vanishes* Hence the 
only isotropic tensor of order 2 is a scalar multiple of S^, 
If u a . m is an isotropic tensor of the third order, 

u 'ikm = (S*j — 6 ,i) (S* t — C fc |) (S wn “ C mn ) u iln (1), 
and therefore, for all values of i y k t m, 

Cii u lkfn + Vki u iim + Vmi u iki = 0 {^)* 

Take % = k = I* Then 


4 C 12 U l‘>m "I" 

4 + (Wi^m 4" = 0 (3)* 

Now put m = 2 so that c mt = 0 * Then 


& 2 I 2 4 % 2 S = 

+ u m ~ 
U m = 0 


(4). 


From the last equation, and by symmetry, u fkfK — 0 if two 
of i , k t m are equal and the third unequal* Then by the 
first, u ikm is also zero if all of i, k t m are equal; and the 
second shows that 

u ikm ~ ~ u kim* 

If in (3) we put m — 1 , every term vanishes, so that ( 3 ) 
holds* 

Now in (2) if i, m are all different, u j):m is zero unless 
j = i , and then % = Hence ( 2 ) holds* It follows that the 
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only isotropic tensors of order 3 are scalar multiples 

of €jjtm- t 

If Ma-mp is an isotropic tensor of order 4, we nave, 
similarly, 

+ CkfW'tjmp + c mi u ikip + C Pi u ikmS “ ® 

There are only three possible values for i, k , m, p, and 
therefore at least two of them must be equal* We may con¬ 
sider separately the eases where {a) two are equal and the 
other two unequal, (b) three equal, (c) two equal and the 
other two equal, (c£) all four equal. 

In case (a), take i = k = l t m = % p = 3. Then 

C 12 1%23 + 3 + + C 13 'W 132 3 

+ c^« 1113 + + C&Umz = 0 ( 2 )* 

Hence, by the antisymmetrical property of c ikf 

t%2» + ^223 ” %m = 

^3123 + ^1323 — = 

%133 - %122 = 0 

Other instances of case (a) can be obtained by inter¬ 
changing suffixes that are not already equal, and by turning 
the axes so as to bring 3 into the position of 1, 1 into that 
of 2, and 2 into that of 3* Thus (4) gives 

^1133 = ^31122 = Unm ~ u mi = %*22 = ^3311 (^)“ 

And also 

Uj 313 = U 12 12 = 3*2323 = ^2121 — U 3232 = % 131 (®)» 

^3113 1=1 U 2 in ™ %2£3 = U mi “ ?( 2332 = %331 (')* 

In case (6), take i = k = m = 1, p — 2. 

^12^2112 + ^03^3112 "h C 12^l212 + C 13^131t 

+ %%122 + £ 13^1132 + C 21 ^ 11 U + C 23 % U 3 = ® t ^) 4 

The last term shows that 

% ns = 


(3) , 

(4) , 
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and therefore, by interchange of suffixes, all components of 
class (6) are zero* Also, from the coefficient of 

^sna + %aia + %iaa = ^ (10)* 

But in (3) the last term vanishes and we infer 

%I23 "I" ^1223 = 0 (1 f)> 

and therefore 

«1313 + «3112 = 0 (12) s 

whence, by (10), 

% 133 =0 (13). 

Thus all components of class (a) are also zero. 

The coefficient of c 12 in (8) gives 

^ul = ^2112 + ^1212 "J" ^1122 (l^)r 

so that the components of class (rf) are expressible in terms 
of the three types of class (c). 

No further information is got by transforming com¬ 
ponents of classes (c) and (d). Tims if i = k = 1, m = p = 2, 
replacing i or k by j will give a zero component unless j is 
equal to 1; and then the factor c if or c k4 is zero, and the 
relation holds automatically. Similar considerations apply 
if all of i, k , m t p are equal. 

We may denote the components of type (5) by A, those 
of type (6) by p + v t and those of type (7} by ju - v. Then 
(14) gives 

^Ittl = ^$222 = ^3333 = A + 2/4 (15). 

There appear therefore to be three independent isotropic 
tensors of order 4, obtained by taking each of A, /*, v in 
turn equal to 1 and the others to zero. 

In the A tensor, u ikmv — 1 if i = k and m = p, and in all 
other cases is zero. It is therefore equivalent to S iJt 8 maf , 
which is obviously a tensor of order 4, being the product of 
two tensors of order 2, 

In the /4 tensor, u ikmv = 1 if i = ?n t Jc — p t or if i = p, 
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k = m t and % ^ k. If also i = k r the component is 2. Other 


components are zero* This can be written 

^ihmv = "h (^)j 

and is obviously a tensor of order 4* 

In the v tensor, u ikmv = I if i = m, k = p t and = — 1 if 
i = p, k = m, and in all other eases is zero* If also i = k } 
u ikmp is zero. In this case, therefore, 

u ikmp ~ $im B kP — Sfcm (1 J ) * 

This can also be written 

u ikmv = *ijk € mip 0 ®)* 

for if i = 1 ,k= 3, e m = 0 unless j = 2 and then = 1. But 
then € miv = 1 if m= 1, p = 3, — 1 if m= 3, p= 1, and 
otherwise = 0* Thus 

^1313 “ I > %3Sl = “ I 


with corresponding values for the other components* 
Evidently {17) and (18) represent a tensor of order 4* It 
has already appeared in Chapters i and vi. 

The general isotropic tensor of order 4 is therefore 

H- p (SjfnSjtj, + (-0), 

where A, /w, v are scalars* 


EXAMPLE 

SfJfeS m p10fjfc 131 & m pWiit 


Prove that 


CHAPTER VIII 


ELASTICITY 


In an clastic solid, as in a fluid, the distance between any 
two particles of the body usually varies with the time. 
The body, however, has an equilibrium configuration that 
could persist if the external forces were zero or constant. 
We may take this as a standard of reference* If a 
particle actually at x t would be at x t — u t in the standard 
configuration, we call u { the displacement at x { ; in practice 
the squares of the u t can usually bo neglected* Evidently 
% is a vector* 

At a point x t + y it where is small, the displacement is 


%i + v it where 


du f 
dx k ' 

~ ( e ik £ik) lIk 


V * = 3^ y * 


( 1 ) 

( 2 ), 


where e /k and $ ik are the symmetrical and antisymmetrieal 
tensors 


—*(S+fe)= <3 >- 


If e ik is zero at x it the displacement has the same form 
as that due to a general displacement Uj together with a 
rotation expressed by the tensor £ ik . Also, if e ik is every¬ 
where zero. 


= t JL 

dx m 2 cx m Vca;,. exj 

= i A /A* + + A 

2 dXj \3;t m ' exj 2 dx k \3ar„ ^ oxj 

_ _ q 

ax f dx k 


(*)• 
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Hence the are constants and the rotation is the same 
everywhere. The vanishing of the e ik is therefore the neces¬ 
sary and sufficient condition for a pure rotation. 

Now consider the part of v* due to the e ik . If we consider 
the quadric surface 

= r* (5). 

where r is a constant chosen so as to make the surface pass 
through y it the normal at y t has direction cosines pro¬ 
portional to e ik y k and therefore to v ( . The displacement 
due to e ik is therefore parallel to the normal at y f to this 
quadric. This coincides with the direction of the radius 
vector to y, if y t is on a prmcipal axis. There are therefore 
three directions such that the relative displacement due 
to the e,* is in the direction of y it and these directions are 
mutually perpendicular. If we take new axes £ 1 , £2 > £3 
along them, the quadric reduces to 

~ r * w> 

and all terms e,i with j # l are zero. The displacement in 
the $1 direction is now , so that all distances in this 
direction are increased in the ratio 1 to 1 + e^'. The dis¬ 
placement due to the e tk is then the resultant of three 
homogeneous strains parallel to three orthogonal axes. 

We see therefore that the displacement in any small 
neighbourhood can be represented as the combination of 
a rotation with three extensions at right angles. The latter 
express the changes of size and shape of an element of the 
solid. For this reason t ik is called the strain tensor and £ ik 
the rotation tensor. Evidently e ik has six independent 
components. For we can have 

= ey t , v z = v 3 = 0 (7), 

making all the e ik zero except e n = e; similarly e** and % 
may exist independently of the others. Also, if 

v x = 0 , v t = ez, v 3 = etj 


( 3 ), 
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all the e ik vanish except = e S2 = e. Similarly % and e ia 
can be assigned independently. 

In an elastic solid the internal force across an element of 
surface is in general inclined to the surface. If the area of 
the element is dS, the force across it must be specified by 
three components of the form p, u dS parallel to the axes; 
where n may be regarded as indicating the normal to the 
surface. If we consider a small parallelepiped with edges 
dx t , dx a , dx 3 , centred at (s t , x 2 , x 3 ), the force across the face 
of area dx 2 dx 3 centred at x x + {dx 1 is (p u , p lt} p ls ) dx t dx 3 , 
where the p u are evaluated at x t + ldx 1 . The force on the 
opposite face is — {p tl , p n , p n )dx z dx 3 evaluated at 

x \ — and the total is c ^ i dx 1 dx 2 dx 2 * In general the 

force in the x t direction due to the stress across the faces 

of x k constant is ^ - dx 1 dx i dx at and we take account of 
cx k 

all faces by using the summation convention. If the 
acceleration of the element has components /, and the 
density is p r the mass is pdx x dx 2 dx ^; while if the bodily 
force acting has components X f per unit mass* the equa¬ 
tions of motion are 

° ~dz‘ + px * < 9 >* 

Ihe system of quantities p ih constitutes a symmetrical 
tensor of the second order. To prove this, consider first a 
plane whose normal has direction cosines intersecting 
lines through parallel to the coordinate axes at short 
distances from thus a small tetrahedron is formed, 
whose sides are of order l , say. Let the area of the sloping 
face be dS; then those of the others are a^dS. Consider 
now the rate of change of momentum of the matter within 
this tetrahedron. Evidently the contributions from the ac¬ 
celeration and Xt are of the order of the volume, that is, 
of R The force across dS has magnitude That across 
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the face of x k constant is — Pm times the area of the 
face, that is, - p k( a ki dS. But dS is of order IK Hence we 

‘“’ ,e (io), 

and hence if l is indefinitely small we have for the stress 
across a plane normal to at x if 

Pu “ a uPki ( n )- 

Now consider three perpendicular directions with direction 
cosines a, t (j = 1, 2, 3). The force per unit area across a 
plane perpendicular to one of these axes, in the direction 
of x { , is given by pn . Resolving this along the direction of 
one of the new axes x t we have, therefore, 

Pil = <hi a ilPki = a U a *lPik 

by interchanging i and 4; which is precisely the law of 
transformation of a second order tensor. _ 

Consider again a small parallelepiped centred at Xi , with 
edges parallel to the coordinate axes, and form the equa¬ 
tion of moments about its centre. The contributions from 
/* and X* are of order l* at most, where the edges have 
lengths of order l. The moment about an axis parallel to 
of Hie stress p ai in the face x % = constant and parallel to *, 
is the product of pa into the area of the face and the dis¬ 
tance of the face from the centre; that is, to order l 3 , 
lp a dx 1 dse t dx t . The opposite face makes an equal contribu¬ 
tion. The stress p 12 in the face x t = constant and parallel 
to x k tends to turn in the opposite direction. We have, 

therefore, ^ __ ^ dx 1 dx 2 dx 3 = O (i 1 ) (13), 

and therefore when we make l indefinitely small wo must 

h&VG / | A V 

Pn = Pn < 14 >’ 

and in general 

Ptk = Pki 

so that the tensor p ik is symmetrical. 
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Now consider the energy interchange between the small 
parallelepiped and its surroundings. The stresses across the 
face of area dx 2 dz z centred at | dx l are doing work on 

the element at a rate dx^dz^ and the contribution 

from the two opposite faces is 

9 

6 - (p kt n k )dx 1 dx. i dx< t . 

Thus in all the stresses are doing work at a rate 
dx k 

The external forces are doing work at a rate pX*ti f dr* The 
kinetic energy of the element is ^pu/dr, and is increasing 
at a rate pUjfidr, (We consider the actual specimen of 
matter occupying the element of volume dr at time t ; thus 
its mass is pdr , and if we keep to the same piece of matter 
at time t + dt the mass is unaltered* If we considered the 
variation of energy within a given element of volume we 
should have to allow for the variation of p and the fact that 
the matter moving out of the element is taking its energy 
with it.) The rate of performance of work on the element 
therefore exceeds the rate of increase of kinetic energy by 

jJ^ (P*(i‘,) + pXt'h - pujir dr 

- {£ <?..*.) + f-M—i, + Oi’-)} 

by the equations of motion, 

= { P “tx) dT (I6) ' 

This work is stored up as internal energy of the element of 
volume* Evidently from its form it is a scalar. 

In any elastic solid the internal energy is a definite 
function of the state of the solid* In any change of state 
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the increase of internal energy therefore depends only on 
the initial and final states and not on the method of passage 
from one to the other. Now we have seen that six of the 
e,* are independent, and if the element acquires displace¬ 
ments 8u ( in time St the corresponding increase of internal 


energy is p ki S dr. There is an apparent asymmetry 

according as i and k are equal or unequal, lhus Pnd T has 
coefficient 8 (dujdxj = 8e u , but t = 1, k = 2 contributes 
p u S (dujdx?), and i = 2, k = 1 contributes p l2 8 (Su^dx,), 
the two together giving 2p w Se ia . But this is the same as 
p n Se lt + p n §e 21 , and the whole contribution from the 
changes of strain is p ik $e ik dT. Also during the process an 
amount of heat SQdr may be absorbed. If then Edr is the 
internal energy of the element. 


82? = p ik Se ( * + SQ {17). 


Since E is a definite function of the state of the system, 
and six of the 8e (i are independent and determine the 
other three, 62? depends on the changes of the temperature 
and of the six independent e ik and has a definite value in 
whatever order these changes come about. But 



/32?\ 

Ur) 

\ce n /zQ^Q 





(18). 


If the absolute temperature is 0, and a certain amount of 
heat 8 Q is absorbed without change of any linear dimen¬ 
sion, the rise of temperature is related to S<2 by the rule 


8Q = pcS0 (19). 

where c is the specific heat at constant strain. If there is 
also a change of strain, since 8 Q and 80 are scalars, we must 
have 

SQ = + />c 80 (-0), 

where the q ik constitute a tensor of the second order. But 8 A’ 
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and 8 Q‘,0 are perfect differentials. Hence if wo replace the 
six independent e ik by e r we can write 


82? = (p r + q r ) Se r + pcSB 

t-iK + 's S9 


( 21 ), 

( 22 ), 


and 


jjj^ (P* + = de t + q ' > ’ cO ^ Pr + ~ ?p ^ 


3e r ' 


1. /fc\ _ A 

de r \6) de,\9)’ 


3 

30 


»K(f) < 2 *>* 


It follows at once that if 0 is kept constant Ep r Se r and 
Y.q r he r are perfect differentials. Also 


3 

30 


(Pr + ?r) = 0 


3e r 





<h 

0 


(25), 


and tlicrcfore 



If 80 = 0 we can write 

Zp T U r = 8 ^ ( 27 ), 

where 

2IK *■= Cq -t- 2c T e, 4- C„e f e. + 0 («») (28). 

The c 0 ,c T) c rs may be functions of 0. Then if we retain 
only terms in IK up to order e a , 

p r = c r + c„e, (29). 

The c r represent the stresses that would remain if the 
strains e T were removed without change of temperature. 
In most practical eases the original state is one of uniform 
temperature and no stress, so that c r =0. If there is a rise 
of temperature 6' under no stress, an element will acquire 
displacements 

V t = a ik 6'y k (30), 
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where a ik is a second order tensor expressing the thermal 
expansion. Thus 

— i i a ik + a*,) 9* — Put 9' {31), 

where fi ik is a symmetrical tensor; and 

Pit = C ik + C ik t mvPmp9* (32), 

where c ikmp is a fourth order tensor. 

But by hypothesis this deformation takes place under no 
stress and therefore p ik =» 0, This determines c* fc , and our 
formula for the stress is 

Pik *= (emp - p mp p) (33). 

The coefficient c ikmp is the coefficient of in IF, 

Since there are six independent e r there are twenty-one 
possible terms in a quadratic form c„e T e 9i and therefore 
there are twenty-one coefficients c (k mf) , They clearly form 
a tensor of order 4; such a tensor in general would have 
eighty-one components, but this satisfies the symmetry 
relations that it is unaltered if we interchange i and £, or 
m and p , or % and k together with m and p together, 

From (26) and (33), 

3 3 

= 9c mp (c^ ffl p) -}- 9 (Gik'tnpPmp $ ) 

3 

“ ‘ 4- fitnv (■* 0 

if 9* is small. The second term does not involve the e ik ; 
the first is small of the first order in the e ik . 

Many solids are isotropic; that is, they have the same 
properties in all directions. This applies to vitreous (glassy) 
solids and to mixtures of crystals oriented at random. In 
that case a uniform rise of temperature in an element 
gives an equal expansion in all directions and 

v t « a9% 


(35) 
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simply; then a is the coefficient of linear expansion and 
Pik 33 (36). 

The second order terms in W% constitute a scalar; and 
we have 


3 Pit 

3<w 


C ik t mp 


(37), 


a tensor of order 4, If it is isotropic it must be of the 
form (20) of Chapter vh. Tlien the linear terms in p {k give 

Pik ~ € (k.mp e mp (33) 

= A3 JJfc 8 mi ,e mjJ 4- /x {S im S* p 4 - & tp 8 m ) e mp 

+ v € mp 

= 4- p {S im e mk 4- Sfp 6 fc p) 4- v (S rf7 , ~ &ip&ki>) 

= AS lfc e mm 4- p (e ik 4- e**) 4- v (e ik — e ki ) 

= + 2/xe fJt (39), 

the last term vanishing since e ik is symmetrical. 

Then 


( 40 ) 

— Ae^e mm + %i ie ik e tk 

= A A s 4- 2 /a( e^H- €fia a H“^ 33 a + 2^“4* 2 e 3l - 4 - 2 c J 2 m } (41), 
where A = e a = du t jdXi (42). 

The scalars A and y represent properties of the material. 
Both are positive, 

YV r e can also write 


2W t = (A 4- 2/a) A 2 4- 4/a (e^ 2 h ej, 2 4- 

“ “ ^U%a) (43). 

This appears to differ from the form in Lovers Elasticity t 
1906, p, 100 , but the present e ik differ from Love's strain 
components. My e n is the same as his, namely 3ujdx l ; but 
my £*23 is only half his, so that his assemblage of strain 
components is not a tensor. 
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If all the e it were equal to one another and therefore to 
JA, we should have 

2W % = (A 4- Im) A 1 (44), 

Pik = (^ + I/*) A (* = j Ptk = 0 (i t £ k) (45), 

In general we write 

A + i/x = k (46), 

and call k the bulk-modulus, 

2IF 2 = (A+f/*) A 2 

+ 2/4 (fin* + && 2 4- 4 4 2ey 2 ) 

~~ Im (fin 2 + 4- figg" 4 + ^33 e u 4 

-(A 4 1/4) A 2 

+ (^lt 2 4 e 23 2 4 % a ^ e 22^33 “ ^33 fill — C 11 

4 4/X (fc^ 2 4 Sin 2 + e ! 2 2 ) 

- (A 4 I/ 4 ) A 2 

4 5/4 {{#22 “ %)“ 4 (% fin ) a 4 (fin “ fiia) 2 

4 6^4 6^*4 6^} (47), 

The coefficient of /x vanishes if and only if the strain is 
a symmetrical expansion, and may therefore be called the 
distortional strain energy. 

If we allow for variations in temperature, 

Pik ™ A (A — 3a8 ) S fi 4 2/4 (e^ — 

= {AA - (3A 4 2 m ) aB*} S ik 4 2pe ik (48) 

and ( 49 )- 

Every term in p ik is of the first order in the displacements; 
but cB'/dd l and therefore gives rise to a constant term. 
This term is 


— 0 * *ikcc 


(60)- 
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If 3(3 = 0, so that no heat is lost or gained by conduction, 
pcSd = — qik$ e ik — ~ 3fca08 iJt Se (fc 

= - 3fco0SA (61), 

and therefore, if the strain takes place adiabatically, 


and 


Ptk 




e-_Maf a 

P C 

9k i a i 6\ 


A + —-— j S«A 4 - 2 p,e it 


(52) , 

(53) . 


Thus in an adiabatic disturbance the constant A is in¬ 
creased above its value for a standard disturbance to A', 


where 


A' = A + 


9&»g a e 

pc 


(64), 


while p is unaltered. The bulk-modulus k is therefore also 
increased to k', where 

k — k +- (55). 

pc 

In a simple thermal expansion at zero stress the absorp¬ 
tion of heat S<2 is equal to pc„S0, where c„ is called the 
specific heat at zero stress, and is the specific heat mea¬ 
sured in ordinary experiments. Then 

P<yS0 = &Q = pc BO + q ik Se (Jt 
1 = pcBO + 

= pcS0 + 3/.' <x8 .3a 50. 


Thus 


, = c( 


1 + 


9&a 2 0\ _ k-c 
pc / k 


(56), 


The equations of motion at constant temperature, ii the 
properties A and p are uniform, can be written 

pfi = g“ (A 8 a A 4 2p.c fk ) 4 pX t (57) 

3 / d u k duA 


9 , \ * i 0 /0ttj . vUA yt 

'dx t ^ ^ + ^ 5ai* \cx t + 3a: J p * 

■- {(A 4- p) A} + + pX t 


(58). 

jCT 


7 
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If the changes are adiabatic X must be replaced by 

X\ 

If there is any heat conduction , the absorption of heat 
per unit time by the element of volume dr is 



where K is the thermal conductivity. Then the equation 
of heat conduction is 

*%+*»* s?-*™ <59) - 

if we write 

- aS tk 6 f + «(/ m> 

so that e ik is the strain due to the stresses, 

+ (61), 

ct ct ct 

and the equation becomes 

P c, || + Ma6 = K^-6 (62). 

The element dr originally had volume 

0 IT] , % 2 ^2 * 

d (%, x 2 , ^ 3 ) 

-{ i -(S + 0 S + S) +O( '“ > i <ir<63) ’ 

so that its density was p {I + A). Hence 

dA = _ 1 dp 
dt p dt 

This is the equation of continuity. 


(M). 


CHAPTER IX 


HYDRODYNAMICS 

In comparison with a typical elastic solid, a real fluid shows 
a great resemblance and a fundamental difference. The 
force per unit area across an element of surface parallel to 
a coordinate plane again constitutes a symmetrical tensor 
of order 2, and for the same reasons. The equations of 
motion are still (9) of Chapter vm, and also the rate of 
performance of work on an element of volume dr still has 
the form (16). The difference is that the internal energy in 
a fluid does not depend directly on how much it has been 
deformed. However the fluid is moved about and stirred 
up, provided it returns to its initial position, density, and 
temperature, the initial and final internal energies are 
equal. The deformation, however great it may be, makes 
no contribution; the stresses do work on each element, and 
thereby supply energy, but this is removed in restoring the 
original temperature. If energy of deformation existed in 
a fluid, all particles of it would have a tendency to return 
spontaneously to their original positions when stresses are 
removed, and they have none. Accordingly, while the rate 
of performance of work on the element of volume dr in time 
dt is still p ki du i jdx k dr J where the symbols have the same 
meanings as in elasticity, wo can no longer assert from this the 
equation (18) of Chapter viu, because the change of internal 
energy is not determinate when the changes of the e tk are 
given. The e ik may be as great as we like, but the energy 
does not increase indefinitely apart from changes in density 
and temperature; and the fluid moves in the same way under 
the same external forces whatever its previous history. We 
may say that an elastic solid has a memory; a fluid has none. 
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The stresses in a fluid are related, not to the total de¬ 
formation, but to the rate of increase of deformation. In 
the former notation these have components 

(dn t . 3'/* 


de ik 
J* or 
ct 


1 

$ \dx k dxj- 


The velocities now appear, instead of the displacements 
from an initial configuration, and we now denote the 
velocities, instead of the displacements, by 
We also write 

«--*©-£) (i »' 

so that e ik now denotes the rate of increase of strain and 
£ ik the local angular velocity, usually called the vorticity. 
Now we say that p ik is linearly related to the new e ik ; and 
therefore 

Pik ~ G ik,inv&mv 

where P ik is a symmetrical tensor of order 2 that possibly 
does not vanish with the e ik , and c ik mp is a tensor of order 4, 
Further, in a fluid at rest the stress is an isotropic tensor. 
Hence P ik is isotropic and must be of the form — ph ikt 
where p is a scalar, and c ik mp must be of the form (20) 
of Chapter vn. On carrying out the summation with regard 
to m and p the term in v disappears and we have 

Pik =■ - P 5 ik H~ (®)» 

where A and p are scalars, at present undetermined. 
Evidently p and A are not both necessary to preserve the 
form, and we may introduce the further convention 

Pa - - W» 

so that - p is the mean of the three p ik with equal suffixes. 
This gives 

— 3p =3 — 3p + 3Ae mm + (a), 

and therefore A — — §/* (®)* 

If we consider the internal energy as a function of the 
density p and temperature $ alone, we may consider the 



HYDRODYNAMICS 85 

energy change in a symmetrical expansion. If the density 
increases by S p, there is a contraction in all dimensions in 
the ratio JS p/p and the stresses do work 

(7), 


i Sp Sp 


At the same time there may an absorption or generation 
of heat; then the energy change is 

( 8 ). 


SE = p 8p + 8Q 
P 


and the heat absorbed may be written 

S Q - MSp + pcSe (9), 

where M is unknown and c is the specific heat at constant 
volume. Then 

S E= + MJ 8p + petti (10), 


SQ M « I sp 

e = jSp + s w 


(ii). 




(13). 


and the condition that these quantities are perfect differen¬ 
tials gives 

*--»»© < 12 >' 

i (?) 

80* V 

Tt appears also that p must be a function of p and 9 
alone, for a given material; it does not involve dpjdt>* YV e 
may call p the pressure* 

* This amounts to saying that there is no dissipation of energy in a 
symmetrical compression or expansion. This is true in a gas on the older 
kinetic theory; but Enskog has shown {Kungh Svziutba A hod, IlandUngar, 
no, 4, 1922, p. IS) that p can with greater accuracy be given by 

p ~Jl p9 + idi' 

where ij is a “second " coefficient of viscosity. But yjffx, is only of the order 
of tho square of the ratio of the volume of the molecules themselves to the 
whole volume of tho gas. 

In a liquid the coefficient of dp/dl t if it exists, is within the experimental 
error; ao analogous statement is true for imperfectly elastic solids. 
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The stress components may now be written 

Pth = — (?>+ Si* + (14), 

where p is a function of the density and temperature and 
yt. expresses a property of the fluid. We call p the coefficient 
of viscosity. For a uniform fluid we have the equations of 
motion 

& = dX k {“ + fr 6 ™) S« + fl + g^)} 

= - g|. {p + fjUC fflm ) + fiWUi + IX ~ e kk 


" ” fa f (? ~ 4^) + 

whore we now write the scalar 

e m - = A 


(15), 


( 16 ). 


In time dt the outflow from a volume element dx x dx 2 dx^ 
is {3 {pu,)jdx ( } drdt. The mass within the element therefore 
decreases at a rate {3 (pu ( )/dx ( } dr, and we have the equation 
of continuity 

dt = ” ete ( {pUi) (17> ’ 


or 


da Bn* jl 

dt ^ 3a:, ~ ~ P A 


( 18 ). 


Here ^ denotes differentiation with regard to the time 

following a given particle of the fluid* so that ^ are 
given by 


and 


dxj 

—5 = Uj 

dt 1 

d _ d dx t 3 _ 3 3 

di = dt + St dx t ~ 0 t +Ui dx i 


(19), 


( 20 ). 


In most works on hydrodynamics this operator is denoted 
by DjDt, but I see no reason for departing from the usual 


HYDRODYNAMICS 


87 


notation for total differentials, since this particular one 
is the only total differential that occurs in hydrodynamics. 
The acceleration components are given similarly by 



3 u t „ du t 
Ti + Uk Bx k 


( 21 ). 


Consider now the circulation Q about a closed contour* 
We have 


w 1 ■ & \ 0 uidxi= \M dxi + ic u,dut 

" \c ( x<+ \ a*f) dx<+ \c UidUt (22) ‘ 

The integral [ u { du { always vanishes; X t dx ( vanishes 

Jc J C 

for bodily forces derived from a potential Then 

-£ = \dk { - 

= ~\cl dx { {P + M dXi + \cl <k (2lie ‘ k) dXt 

= - jj d (p + + jj ^ {** (- 2f« + 2 §J)} 

" - f 1 d (p + §/*A) - { - Jj~ dx t 
Jc P Jc P ax k 

+ [ - If <*«*+ f (23)* 

Jcp Q%k JoP 

If /X = 0 this reduces to the circulation theorem of classical 
hydrodynamics. In many problems of real fluids p. is 
small and constant, and A small. Then the integrals 

f - d (§/iA), [ - f - du k , and f - pdA are zero or products 

JcP Jcp jgp 

of two small quantities, and can be ignored. - dp is zero 
if p is a function of p, which is true in many cases, though 
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exceptions arise when the temperature or the composition 
varies from place to place. In the commonest case, to 
considerable accuracy, 

f— lk { ^‘ )dXl !24) ' 

and circulation arises only from variations of vorticity in 
the neighbourhood of the contour. It follows that in a 
fluid originally at rest or in irrotational motion circulation 
can arise only through the diffusion of vorticity inwards 
from a boundary. 

The work done on the element dr per unit time exceeds 


the rate of increase of kinetic energy by 



du { , , 

Pkt dr — Pik&ikdr 

(25), 

and 

Pik&tk — — {P + fp-A) S<* 6 <k + e ik & ik 



“ — (P + f j^A) A + 2 {ie tk e tk 

(26). 

But 


(27), 


and if p is a function of p (including a constant as a par¬ 
ticular case) this is a differential with regard to the time 
and yields on the whole zero if the original density is ever 
recovered. This term may therefore be considered to give 
the increase of internal energy due to compression. The 
remainder may be written 

<5 = - + 2pe ik e ik 

«= I** {(^a - «*3) 2 + {% - e n ) s + fat - %,)*} 

+ 4^ (e n * + + <h 3 ») (2 $). 

Thus fi> is analogous in form to the distortiona! strain 
energy of elasticity ((47) of Chapter vm). It is essentially 
positive, and therefore represents work doneon the fluid and 
continually stored as internal energy. In (8), therefore, if the 
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change is an actual one, the term (27) contributes thephpjp, 
while (28) contributes to SQ. So long as the initial and final 
states are given, it is immaterial whether SQ represents 
heat conducted into the element, absorbed from radiation, 
or generated chemically within it, or mechanical energy 
dissipated into heat by viseosity. 

We notice that <l> can vanish if , and only if, 

= &22 = 

% = 631 = £12 = 0 ( 29 ), 

so that the deformation represents a symmetrical expres¬ 
sion or contraction. 

If we consider any finite volume, the rate of dissipation 
within it is 

JJJ cPtfr = jjj (- 1 /iA* + 2fj.e tk e ilc ) dr (30). 

But e (k e ik = UU+llf£ k (31), 

and 



(32), 


where l k is a direction cosine of the outward normal to the 
boundary. But 

= \\n k u i t ik ds + i\\ dl J ) ?as 

- J f U tkm uM m dS + 1 dS (33), 
where djdn denotes differentiation along the outward 
















90 


HYDRODYNAMICS 


normal, q is the resultant velocity, and is the vorticity 
vector associated with Also 


( 34 ), 


and in all 



ft- j" [ dS + 4/* || 


+ V|}j(- JA>+f„f „ + §-*)* < 35 >' 

This form is useful in such a problem as that of waves 
on deep water, where the viscosity is small; if it were 
absent we should have a permanent oscillation in a normal 
mode. The vorticity is negligible except within a distance 
from the bottom of order (vjy)$, where fx = vp and 2?r/y is 
the period of the motion. Then in the second integral 
is zero at the free surface and Ui is zero at the bottom, so 
that this integral vanishes. A is negligible everywhere. 
Near the bottom % (k contains a factor proportional to the 
velocity that would exist there on the classical theory, and 
thus can be made indefinitely small for deep water. Hence 
the important term in (35) is 


and even this vanishes at the bottom, so that it need only 
be estimated at the free surface. 

If a portion of the fluid is compressed without change of 
temperature, we define the bulk-modulus k by 


1 dp _ 1 
p dp k 


( 36 ), 


and if it expands under change of temperature without 


HYDRODYNAMICS 


91 


change of pressure we define a coefficient of volume 
expansion a by 

1 

P 

Then for small changes of temperature and pressure 

S P - P (~£ - «5fl) (38). 

Also in a free expansion under constant pressure 

SQ = MSp + pc89 

= {« | (?). ap + pc | m (39), 

where the partial differentiation is to be carried out at 
constant density. But this makes 

^ - olIc ( 40 ), 


30 


= — a 


(37). 


and therefore 



a 2 * 01 30 


(41). 


so that c v , the specific heat at constant pressure, is given by 


c P = c + 


a*k$ 


P 


(42). 


This form is analogous to that of (56) of Chapter vm, the 
present a being the coefficient of volume expansion, and 
the previous one that of linear expansion. 

In an adiabatic change SQ = 0; then 


and therefore 


so that 


-ak$^ + pcse=0 

p 

ad Bp = {pc + a 2 k0) 80 


*-*( 1 + ir)? 


( 43 ), 


( 44 ). 
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Thus the bulk-modulus for adiabatic changes is 

F=fe( 1 + ^) = ^ ? (45). 

\ pc / c 

The equation of heat conduction needs to be modified 
to allow for heat generated Internally, In time §£* per 
volume dr 7 we have 

+ £;(*£) 81 

— pcSd 4* ak8A$t (4G) f 

whence pc jj| 4- ct&flA = ^ gjj ( 47 )- 


EXAMPLES 

1, Obtain the equations of motion in terms of the stress 
components by considering the momentum of a finite volume of 
any form and applying Green's theorem, 

2* Similarly obtain equation (25). 


I N D E X 


Acceleration, 24 
Angular momentum, 35 
Area, plane, as vector, 20 
Axes, change of, 1 
dynamical, 24 
rotating, 40 
Axis, central, 45 

Centroid, 33 
Change of axes, 1 
Circulation, 50, 87 
Couple, 43 
Carl, 14, 53 

D*Alembert's principle, 33 
Dielectrics, 52 
Direction cosines, l, 16 
Dissipation, 38 
Divergence, S, 53 

Elasticity, 71 
Electromagnetism, 61 

Electrostatics, 50 
Energy, kinetic, 25 
potential, 26 

Equation of continuity, 58, G2 
Force, 25 

line of action of, 42 

Gauss's theorem, 49 
Gravitation, 49 
Greeks theorem, 48 

HumHton'a principle, 23 
Heat conduction, 32 
Hydrodynamics, 83 
Hydrostatics and classical hydro* 
dynamics, 56 

Inertia, momentsand products of, 37 
Isotropic solids, 73 


Lagrange's equations, 23 

Magnetism, 53 
Moment of force, 43 
of vector, 35 

Null planes and points, 46 

Parallel, 17 
Potential, 49 

Rigid bodies, 29 
Rotating axes, 40 
Rotation, 21, 30 

Scalar product, 7 
Scalars, 3 

Specific heat, 81,91 
Stokes's theorem, 48 
Strain tensor, 72 
Stress tensor, 73 
Suffixes, dummy, 5 
Summation convention* 3 

Tensors, 2 

autisymmetrical, 5, 3 
isotropic, 66 
symmetrical, 5 
&ik* 6 
*ikm > ^ 

Vector product, 7, 14 
Vectors, 2 

with given divergence and curl, 
68 

Velocity, £4 
Virtual work, 34 
Viscosity, 86 
Vorticity, 60 

Work function, 26, 50 
Work, virtual, 34 
Wrench, 45 


Kinetic energy, 25 
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